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ABSTRACT 

The observed bimodality in radio luminosity in galaxy clusters is puzzling. We in- 
vestigate the possibility that cosmic-ray (CR) streaming in the intra-cluster medium 
can 'switch off' hadronically induced radio and gamma-ray emission. For self-confined 
CRs, this depends on the source of MHD wave damping: if only non-linear Landau 
damping operates, then CRs stream on the slow Alfvenic timescale, but if turbulent 
wave damping operates, super-Alfvenic streaming is possible. As turbulence increases, 
it promotes outward streaming more than it enables inward turbulent advection. Cu- 
riously, the CR flux is independent of V/ (as long as it is non-zero) and depends only 
on plasma parameters; this enables radio halos with flat inferred CR profiles to turn 
off. We perform ID time-dependent calculations of a radio mini-halo (Perseus) and 
giant radio halo (Coma) and find that both diminish in radio luminosity by an or- 
der of magnitude in several hundred Myr, given plausible estimates for the magnetic 
field in the outskirts of the cluster. Due to the energy dependence of CR streaming, 
spectral curvature develops, and radio halos turn off more slowly at low frequencies 
- properties consistent with observations. Similarly, CR streaming rapidly turns off 
gamma-ray emission at the high-energies probed by Cherenkov telescopes, but not at 
the low energies probed by Fermi. CR mediated wave-heating of the ICM is unaffected, 
as it is dominated by ~GeV CRs which stream Alfvenically. 

Key words: radiation mechanisms: non-thermal, turbulence, galaxies: clusters: gen- 
eral, radio continuum: general. X-rays: general 



1 INTRODUCTION 



Cosmic rays (CRs) in the intra-cluster medium (ICM) can 



arise from structure formation shocks (Miniati et al. 2001 



Pfrommer 



[20081, turbulent reacccloration of existing non- 
thermal particles ([Brunetti fc Lazarian 2007[), galactic winds 
and supernovae (Volk, Aharonian & Breitschwerdt 19961, 



and radio galaxy jets ( Enfilin et al.|1997 Ensslin et al.|1998 



[McNamara fc Nulsen||2007[ ) . They are visible in clusters in 
radio emission, and gamma-ray emission (via hadronic inter- 
actions). However, unlike in our interstellar medium (ISM), 
CRs in the ICM are energetically subdominant; for instance, 
current upper limits on CR-induced gamma-ray emission in 
Perseus suggest CRs are ;$ 1 — 2% of the thermal energy 
density (lAle ksic et al.|[2012 l. Why then are CRs in clusters 
of astrophysical significance? Firstly, unlike in the ISM, cos- 
mic ray protons (CRp) with E < lO'^GeV remain confined 



and have lifetimes of order a Hubble time ( Volk, Aharonian 



fc Breitschwerdt 1996 Berezinsky, Blasi fc Ptuskin 19971; 



they therefore encode archaeological information about the 
cluster assembly history as well as AGN and supernova ac- 
tivity. Secondly, the ICM provides stringent tests of plasma 



physics in a regime very different from the ISM. CRs in 
clusters represent an opportunity to study the unknown ef- 



ficiency of shock acceleration (Blandford fc Eichler 19871 
in a low Mach number TVf ~ 1 — 5 and high plasma beta 
/3 ~ 100 regime. Transient radio phenomena can also teach 
us about magnetic field amplification at shocks. Thirdly, 
even a low level of CRs could have interesting astrophys- 
ical implications. These range from pressure support (thus 
affecting the use of clusters for cosmology) to heating which 
suppresses cooling flo ws ^^Guo fc Oh^^2008) or energizes fil- 



aments (Ferland et al.|2008 20091, and distributing metals 



dran , 



and heat via buoyancy-induced turbulent convection ( Chan- 
Rasera|2007l [Sharma et al.|2009[ |. 



One drawback of cosmological simulations of CRs in 
clusters is that they generally do not include CR transport 
processes; the CRs are assumed to be frozen into the gas, 
and advected with it. In practice, CRs can move relative 
to the gas by streaming along magnetic field lines down a 
CR gradient, as well as diffusing across field lines by scat- 
tering off plasma waves. As CRs stream, their momentum 
anisotropy excites plasma waves, which in turn scatter the 
CRs, isotropizing the CR distribution in the frame of the 
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waves. This generally limits streaming speeds to the speed 
of the waves, which is the Alfven speed da- In our ISM, 
rapid p itch-angle scattering due to the CR streaming insta- 
bilit^f] ([Lerche]l967| p<ulsrud fc Pearce|l969||Wentzel|l969 



Skilling 1971 1 can explain the observed spatial isotropy of 



CRs, as well as the escape time of CRs from the Galaxy 



I Schlickeiser 2002 Kulsrud 20051. Applying the same CR 



self-confinement scenario to the ICM, the low implied drift 
speed of CRs vd ~ va ~ lOOkms"^, seems to justify ne- 
glect of cosmic ray transport. Early calculations of CRs in 



isolated clusters (Boehringcr & Morfill 1988 Loewenstein, 



Zweibel & Begelman 1991) which considered CR diffusion 



at a level comparable to ISM values found it to be negligible 
as a transport process. They argued that if CRs were in- 
jected at the cluster center by an AGN, they would quickly 
dominate pressure support at a level inconsistent with obser- 
vations. Guo & Oh ( 2008 1 resolved this in their calculations 



of CR heating by allowing CRs to be transported by rising 
buoyant bubbles, as seen in high resolution Chandra im- 
ages, which are subsequently shredded by Kelvin-Helmholz 
and Rayleigh- Taylor instabilities to disperse the CRs. 

However, this assumption of slow CR streaming and dif- 
fusion may not be fully justified. The plasma waves which 
scatter the CRs are also subject to a variety of damping 
mechanisms. If damping is stronger in the ICM than in the 
ISM, then pitch-angle scattering of the CRs off the attenu- 
ated waves will be reduced, and the CRs retain some mo- 
mentum anisotropy in the frame of the waves. They can 
therefore stream faster than the waves, and will no longer 
be limited by the Alfven speed. In principle, if the waves are 
very strongly damped, the CRs could stream up at speeds up 
to ~ c. While the possibility of super- Alfvenic or even free 



streaming was appreciated early on ( Kulsrud & Pearce 1969 



Kulsrud fc Cesarsky|1971||Skilling|1971| ), generally Dp ~ va 
and diffusion coefficients appropriate for the ISM have been 
uncritically applied to the ICM environment. In an influen- 



tial recent paper, Enfilin et al. (20111 noted the interesting 



possibility of super-Alfvenic streaming in the ICM, adopted 
the sound speed Cg as a characteristic streaming speed, and 
were the first to discuss the wide-ranging observational con- 
sequences. 

A particularly interesting possibility they focused on 
was whether the interplay between advection and stream- 
ing could be responsible for the observed bimodality seen in 
radio halo luminosity. Giant radio halos are generally only 
seen in disturbed clusters which show signs of merger activ- 
ity. This bimodality has been a stumbling block for hadronic 
models (e.g., Pfrommer fc Enfilin] |2004J ). These models 
track the long-lived CR protons (CRp) formed during struc- 
ture formation shocks, and find that the secondary electrons 
formed when the CRps undergo hadronic interactions are 
sufficient to explain radio halo observations. Transience or 
correlation with turbulence is generally not expected in such 



^ In principle, CRs can also scatter off MHD turbulence, though 
this is thought to be weak due to the increasing anisotropy 
at small scales, with power concentrated in modes with wave- 
vectors transverse to the B-field, while CRs efficiently scatter off 
the parallel component l |Chandran|200()| |Yan fc Lazarian 2002 ) . 
Fast magnetosonic modes could potentially scatter CRs more ef- 
ficiently ( [Brunetti fc Lazarian|p007[ l , but a treatment of this is 
beyond the scope of this paper. 



models. As a result, radio halos are often attributed to the 
re-energization of seed electrons by Fermi II acceleration 
when the ICM turbulence becomes transonic during merg- 
ers ( Brunetti et al.|200T Brunetti fc Lazarian|2007 l. As the 
turbulence dies away, the CR electrons (CRe) cool via syn- 
chrotron and inverse Compton emission on a relatively short 
(~ 10** yr) timescale. However, the origin of the seed elec- 
trons is uncertain; low energy electrons will rapidly Coulomb 



cool in the dense cluster center. Enfilin et al. (20111 sug- 



gested instead that transonic turbulence advects CRs from 
the plentiful reservoir on the cluster outskirts. Hadronic in- 
teractions of the inwardly advected CRp with the dense 
central ICM can then produce CRs[^ Once turbulence dies 
down, subsequent outward CR streaming switches off the 
radio halo. This rapid outward streaming also explains why 
radio halos turn off in the original hadronic scenario, which is 
otherwise difficult to understand. For these explanations to 
work, the CR streaming timescale must be relatively short, 
or ud ^ Va- Enfilin et al. (20111 adopted vo ~ Cs and ex- 
amined its implications, but only justified this assumption 
qualitatively. 

This paper aims to critically examine the possibility 
that super-Alfvenic streaming could play a crucial role in 
CR transport in the ICM, by building more quantitative 
models to clarify its plausibility and importance. It has three 
main goals: (i) a quantitative calculation of CR streaming 
speeds in quasi-linear theory and its dependence on plasma 
parameters, when a variety of wave damping mechanisms are 
at play. In particular, we consider the effects of non-linear 
Landau damping and turbulent damping. We also give ex- 
pressions for parallel diffusivities and CR heating rates in 
this regime. Of particular interest is the countervailing ef- 
fects of turbulence: it affects CR transport both by turbu- 
lent advection of CRs from the cluster outskirts, as well as 
damping of CR generated waves, which enables faster out- 
ward streaming. We assess their relative importance, (ii) A 
reevaluation of CR heating due to central injection by an 
AGN ( Guo fc Oh|2008" l, taking these streaming effects into 
account, (iii) A ID simulation of radio halo turnoff due to 
CR streaming, to establish if the required dimming by at 
least an order of magnitude can take place within a reason- 
able timescale. We also calculate how the gamma-ray lumi- 
nosity evolves with time. A non-linear, time-dependent cal- 
culation is needed since the streaming speed itself depends 
on CR energy density. In their analytic calculations, |Enfilin| 
|et al.| ( [2bll[ ) consider a steady-state profile where inward tur- 
bulent advection and outward CR streaming are in rough 
balance. This scenario seems somewhat unlikely; it seems 
more probable that at a given point in time, either inward 
advection or outward streaming dominates. To calculate ra- 
dio halo turnoff, we consider the latter. These calculations 
also enable us to compute a fundamental prediction of this 
model: spectral steepening and frequency-dependent dim- 
ming which arise from energy-dependent streaming speeds. 
We compare these with observations. 

The outline of this paper is as follows. In ^ we calcu- 
late in quasi-linear theory cosmic-ray streaming speeds when 
different wave damping mechanisms are dominant, and de- 



^ Alternatively, low energy relic CRe advected from the cluster 
outskirts could provide seeds for turbulent reacceleration. 
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rive expressions for the resulting parallel diffusivity, as well 
as turbulent diffusion rates. In ^ we describe the equations 
we solve numerically with ZEUS, focusing in particular on 
the CR transport equation. We describe our initial condi- 
tions for the cosmic ray profile, which are tuned to match 
radio halo observations for the Perseus and Coma clusters. 
We also present a test case of CR heating by a central AGN. 
In ij4] we present the results of simulations of radio halo 
turnoff due to CR streaming. In i|5] we show how aspects 
of these results can be understood analytically. Finally, we 
conclude in ^ 



2 COSMIC RAY STREAMING: 
QUASI-LINEAR THEORY 

In this section, we derive the basic equations we use, in par- 
ticular the streaming speeds and diffusion coefficients which 
are used in the cosmic-ray transport equation. Our treat- 
ment is by design semi-quantitive rather than fully rigorous. 



2.1 Cosmic Ray Streaming 

2.1.1 Resonant Scattering and Wave Growth 

Consider a cosmic ray proton with Lorentz factor 7 prop- 
agating along a magnetic field line of strength Bq with cy- 
clotron frequency fio = eBo/{mpc), gyroradius tl — 'jc/flo, 
and pitch angle cosine /x. Since va <C c, Alfven waves are per- 
ceived by the CR as a spatially varying but time-stationary 
B-field. An Alfven wave is resonant with this cosmic ray if it 
has a wave vector k whose component parallel to the mag- 
netic field fc|| equals the parallel projection of the gyroradius, 
i.e. if the resonance condition 



(1) 



is satisfied. This resonance is a requirement both for the 
wave to scatter the CR and for the CR to excite the wave. 
This condition can be easily understood: if the magnetic field 
changes on a length scale much longer than the projected 
gyroradius, the CR will simply follow the field line adia- 
batically, with no change in pitch angle. If the field varies 
on much smaller scales, the CR will see a rapidly oscillat- 
ing Lorentz force during its orbit and remain unaffected, 
essentially only seeing the background field. At resonance, 
the CR sees a constant field due to the wave, and hence 
a steady force. The fcy portion of the wave is the relevant 
one, since it has a transverse magnetic field SB± which can 
exert a Lorentz force on the iiy component of a streaming 
cosmic-ray. 

If the distribution of cosmic rays in the frame of the 
Alfven waves is completely isotropic, then the effect of a 
cosmic ray traveling along the magnetic field line in one di- 
rection is cancelled by an equivalent cosmic ray traveling in 
the opposite direction, and there is no wave growth. How- 
ever, Kulsrud & Pearce ( 1969 1 showed that even a slight 



anisotropy in the cosmic rays-which naturally arises in the 
presence of sources and sinks-will cause unstable growth in 
the waves, caused by momentum transfer from the CRs to 
the waves in the course of pitch-angle scattering. The result- 



ing wave growth rate is ( Kulsrud fc Pearce|1969 1: 

rca(fc||)~no^^^^^^f^-i) (2) 

In the above, ncB.{> 7) is the number density of cosmic 
ray protons with energies large enough to be resonant with 
the Alfven wave for some pitch angle fi, namely tl > l/fc| 
(though since the CR spectrum falls off rapidly with energy, 
generally fcy ~ 1/'"l), "-i is the ion density in the plasma, 
and vu and va are the cosmic ray streaming and Alfven 
speeds respectively. This expression is derived from balanc- 
ing CR momentum loss with wave momentum gain, but we 
can understand its main features qualitatively. The rate of 
wave growth scales with that for momentum loss for a sin- 
gle CR, p oc pfiici ~ {'ynic?){Q,o/')) oc ilo, i.e. the non- 
relativistic, rather than the relativistic gyro-frequency. This 
has to be multiplied by the fraction of ions which can drive 
wave growth, ncR(> 7)/?^i (non-resonant ions simply pro- 
vide inertia, slowing down wave growth) , and the anisotropy 
which seeds the wave growth (ud/^a — 1). 

The streaming instability causes the waves to grow until 
pitch-angle scattering renders the CR distribution isotropic 
in the frame of the waves, i.e. vy) ~ va- If we assume 
{v-d/va — 1) ~ C(1)) we can estimate the growth time of 
the waves. If we assume that the energy density is CRs is 
~ 10% of the thermal energy density, ecR ~ O.lethcrm, then 
ncK^EcK) ~ 0.1ni{i5i) where (Ecr) ~ GeV (as is true for 
most reasonable power-law momentum distributions-e.g., 
see Fig. 1 of |EnfiUn et~al] ( [20071 )), and (S;) ~ keV are the 
typical energies of CRs and thermal ions respectively, and so 
ncn/rii ~ 10~^. A similar ratio holds in the coronal regions 
of our Galaxy. For ~ fiG fields, Slo = eB/mc ~ 10~^s~^, 
implying from equation ^ a wave growth time of r^l^ ~ 30 
yr, i.e. extremely short. 

The above arguments suggest that self-confinement of 
cosmic-rays is very efficient, and should always reduce the 
streaming velocities vo ~ wa- The general success of the 
self-confinement picture for our Galaxy means that this has 
been uncritically assumed in other environments such as the 
ICM, and/or CR diffusion coefficients scaled to the mea- 
sured Galactic values. In fact, vo, and the associated dif- 
fusion coefficient De depend on the amplitude of the wave 
field 5B/B, which can be calculated by assuming equilib- 
rium between growth and damping. If damping processes 
are sufficiently strong, then 5B/B will be insufficient to ef- 
ficiently confine the CRs, and super-Alfvenic streaming is 
possible. We now examine this possibility. 



2.1.2 Non-linear Landau Damping 

Parallel propagating MHD waves do not suffer any linear 
damping. However, they can undergo non-linear Landau 
damping when two waves A & B of slightly different fre- 
quency interact to form a beat wave. This beat wave can res- 
onantly interact with thermal particles with parallel velocity 
identical to the wave's phase speed, — [uja—^^b) I {kA—kB) 
(for parallel propagating waves, v\\ — va). Particles moving 
more slowly than the beat wave will extract energy from the 
wave (thus damping it), while particles moving faster than 
the wave will add energy to it. For a Maxwellian plasma. 
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typically {df /dv)v^^=VA < 0, and damping dominates The 
high frequency wave gives up energy to a combination of the 
low frequency wave and resonant particles. 

To aid in physical insight, we present a simplified deriva- 
tion of streaming speeds to be expected if non-linear Landau 
damping dominates, before employing the formulae from 



Felice & 


Kulsrud 


20011. 


plasma is 


( Kulsrud 


20051 



The damping rate in a high-/? 



- NL ~ -1/ - 



3B 
l3 



:0.3"- 
/i c 



SB 
l3 



(3) 



How can we qualitatively understand the first relation? The 
wave frequency uj — fcyVA sets the fundamental frequency, 
while interactions involving the beat wave arise to second- 
order in perturbed field strength (SB/B)^. Since the reso- 
nant condition = pw; = va implies that thermal particles 
with /X = VA/vi are resonant, it is clear that the damping 
rate should depend on this ratio. However, the exact de- 
pendence only emerges from a detailed calculation-either 
by calculating the slope {df /dv)v„=vj^, or from the plasma 
dispersion relation ( Foote fc Kulsrud|1979 1. Note that since 
Vi/vA = /3^''^/2, we have Fnl oc P'-''' . In the second relation, 
we use the dispersion relation uj = fc||?;A and the resonance 
condition, equation 



In steady state, the Vlasov equation for CRs is (Kulsrud 
|fc Pearce|1969l ): 



' dz 



d_ 



(4) 



where B = B^z, ~ fl{SB/B)'^ and «z = MC. This ex- 

presses the condition that the net streaming along field lines 
is set by diffusion in pitch-angle. In the limit of strong scat- 
tering, we can expand / = /o + /i + /2 + •■• where fo{p, z, t) 
is isotropic and fi{p, z,t,ii) -C fo, /2 ^ fi- Let us also de- 
fine T = /i//o and the scale height Lz{p, z) — fo/{df /dz). 
Integrating both sides of equation Q with respect to fj, and 
dividing by fo, we obtain: 



If we set 



OF 
djj, 



c 
vL, 



3(UD — VtC) 



(5) 



(6) 



&B_ 



c c ncR(> 7) 2 

VA V\ 3Lz Ui 



1/2 



(8) 



where /o oc p " (note that our result differs from Felice & 
Kulsrud (20011, who explicitly specialize to 7 ~ 5 for the 



ISM from the outset). If we scale to numbers characteristic 
of the ICM, we obtain: 



5B\- 



= 1.6 X 10 



-6 (ng?o)^/^go">/^10^'«-" 
{n'_,y/^B,GTli%Ll% 



(9) 



where T4kcv = (T/4keV), B^,g = (B/lfiG), Lz.ioo = 
(LJlOOkpc), nLg = (ni/lO-^cm-3), n'^^o = n^'^h > 
1)/10~^'' cm~^), 7100 = 7/100, and we have scaled to n = 
4.6. Note that n^^{> 7) = 10"^°7"^ '' cm"^ roughly corre- 
sponds to a CR energy density in equipartition with a ~ /xG 
B-field. The fact that (SB/B)^ ^ 1 self-consistently implies 
that quasi-linear theory is applicable. If we insert this into 
equation ([7|, we obtain for the drift speed: 



va[1 + 0.9 



-7loo-^''M (10) 



Several points should be noted. For these parameters, 
streaming speeds do not significantly exceed the Alfven 
speed for the ~100 GeV cosmic-ray protons which in turn 
produce the 10 GeV CR electrons which in turn pro- 
duce ~GHz radio emission. For Alfven speeds of va ~ 
^ BuGn-T]/-^ , this implies radio halo turnoff times of 

3, which may seem too long. How- 



70kms ^ B^an^ ]/;^ , 
t ~ 1.4GyrLiooB 



1 1/2 



ever, note that Lz,n'^^ will be time-dependent functions 
during the streaming process, so it is necessary to check how 
streaming evolves in a time-dependent calculation. Our re- 



sults should be contrasted with those of Enfjlin et al.|( 2011 1, 
who describe similar estimates based on IFelice fc KulsrudI 
(20011, but do not give explicit expressions. Unlike them. 



we find vu ^ Cs for plasma parameters corresponding to 
observed clusters; nothing in the problem singles out the 
sound speed as a reference speed. Note that all the param- 



eters in equation ( 10 1 are observationally constrained, so 
order of magnitude departures are unlikely. Also note that 
even though Fnl oc /3^''^, there is no explicit /3 dependence 
in ud. 



so that {jj.cT{jj.)) — {vd — ^^a) (i.e., the leading order 
anisotropy in the distribution function yields the net drift 
relative to the frame of the waves), this yields: 



{VB — va) 



3Lz 



55 



A 

3L~, 



(7) 



where A ~ r^iSB/B)"^ is the mean free path. In steady 
state, the wave growth rate (equation ([2|) equals the wave 
damping rate (equation ([3|). Together with equation ([7|, 
this gives us two equations which we can solve for two un- 
knowns, VB and (SB/B)^. The result is: 



^ In a high /3 plasma, {df/dv)^^^—„j^ is relatively flat for electrons, 
while still steep for ions; hence, ions dominate the damping rate 
( Miller|199l| l. 



2.1.3 Turbulent Damping 

Another source of wave damping comes from the highly 
anisotropic nature of MHD turbulence ( [Farmer fc Goldreich 



2004[|Yan fc Lazarian|2002[ ). We adopt the |GoIdreich fc Srid 
har ( [1995 1 theory (hereafter 'GS') for strong, incompressible 
MHD turbulence; an excellent summary can be found in 
Lithwick fc Goldreich (20011. Turbulence in clusters is gen- 



erally incompressible since it is significantly subsonic except 
at the cluster periphery. The strong turbulence regime where 
GS theory is applicable sets in at wavenumbers k <, koM]^^ 
(e.g., see Nazarenko & Schekochihin (2011 1); since Ma 1 



in clusters, the theory is clearly applicable, particularly at 
the small scales relevant for CR scattering. GS theory has 



support both from numerical simulations (Cho & Vishniac 
|2000 Maron fc Goldreich^^2001) and solar wind measure- 



ments (Horbury, Forman fc Oughton 2008 Podesta 2009 
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Wicks et al.|2010 Chen et al. [2011 1. It is anisotropic, with 



VA 



1/3 



All (Ax) 
A, 



ilvIHD 
r \ 1/3 

A, 



(^Ax) 



1/3 



(11) 

(12) 



where Ax is the length scale transverse to the local mean 
magnetic field, Wa^ is the rms velocity fluctuation across Ax, 
A||(Ax) is the length scale parallel to the local mean mag- 
netic field across which the velocity fluctuation is v\j^ ,iMHD 
is the length scale at which turbulence is excited with veloc- 
ity perturbations comparable to the Alfven speed va (i.e., 
with Ma ~ 1), and e ~ d^^/Ax ~ ^'aZ-^'Mhd is the (con- 
stant) energy cascade rate per unit mass. Note that Lmhd 
is defined to be the scale at which Ma = 1; if turbulence 
is already sub-Alfvenic at the outer scale, then it should 



be considered an extrapolation. Equation (111 describes a 



standard Kolmogorov cascade in the transverse direction. 



Equation ( 12 1 indicates that an eddy becomes increasingly 
elongated along the magnetic field. Ay 2> Ax as the cas- 
cade proceeds deep into the inertial range Ax <^ I/mhd. It 
can be derived from the assumption of "critical balance", 
which states that characteristic linear and non-linear inter- 
action times are approximately equal at all scales (e.g., see 
Nazarenko & Schekochihin (2011 1). Thus, the cascade pro- 



ceeds primarily in the transverse direction, and most of the 
power is concentrated in modes with transverse wave vec- 
tors. Intuitively, we can understand this from the fact that 
in MHD turbulence, non-linear interactions arise from col- 
lisions of oppositely directed Alfven wave packets travelling 
along field lines. A wave packet is distorted when it follows 
field lines perturbed by its collision partner; it cascades when 
the field lines along which it is propagating have spread by a 
distance comparable to Ax . Since the magnetic and velocity 
fluctuations associated with Alfven waves are transverse to 
the local mean fleld, the cascade proceeds primarily in the 
transverse direction. 

Turbulence therefore suppresses the waves responsible 
for self-confinement of cosmic rays, since they cascade to 
smaller scales before they have an opportunity to scatter 
CRs. In particular, the small scale transverse components 
injected by the cascade mean that the CR no longer expe- 
riences a time-steady force in its orbit; instead it sees an 
oscillating force which leads to inefficient scattering. For 
these same reasons, MHD turbulence scatters CRs inefB- 
ciently ( |Chandran[2000| |Yan fc Lazarian[2002| . The damp 
ing rate of a wave is simply the eddy turnover rate ( Farmer 
fc GoldreichpOOil ): 

r 

Ax X'l^ 
Growth 



(13) 



where we use equation (111 



rates are highest, 
and damping rates lowest, for the most closely parallel- 
propagating waves, i.e. those with the largest Ax- Even if 
a CR-generated wave starts out as parallel-propagating, the 
turbulent cascade injects transverse components which sub- 
sequently cascade. The amplitude of magnetic field fluctu- 
ations across a scale Ax thus define a minimal aspect ra- 
ticj^(fcx/fc||),nin ~ <5B(Ax)/B ~ vx^/vA ~ (Ax/Lmhd)''" 



This should not be confused with the eddy aspect ratio 



(using equation (111 in the last step). From the resonance 
condition k^^ ~ tl, the smallest possible perpendicu- 
lar wavenumber is fcx.min ~ e^^'^{r-L.VA)~^^^ ■ Inserting the 
largest possible perpendicular wavelength Ax ~ fcx,min into 
equation ( 13 1, the minimal damping r ate for a wave with 
fc|| ~ rl^ is ( jParmer fc Goldreich|2004| : 

rturb.min ~ f-^l ^ ■ (14) 

\ri,VA J 

If in steady state we balance the wave growth rate (equation 
([2|) with this damping rate, we obtain a streaming speed: 



v-D 



VA 



1 + 1.2 



T 1/2 

-'^MHD,100"CR,-10 



n — 3.5Tn2(T^ — 4.6) 

7ioo 10 



(15) 



where Lmhd.ioo = -Lmhd/100 kpc. We also obtain 



{5B/B) 



10" 



At first blush, non-linear Landau damp- 



ing and turbulent damping both seem to give similarly 
slow streaming speeds. However, note that v-a ~ va oc 
('^cr''^' "-cr) f'^r these two sources of damping respectively. 
This difference becomes crucial during non-linear evolution, 
enabling CRs in the turbulent damping case to stream much 
more effectively. 



2.1.4 General Remarks on Cosmic-Ray Streaming 

We have now derived streaming speeds for two different 
damping mechanisms, which depend both on CR energy 
(iiD oc 7"'*, 7^'^ for non-linear Landau damping and turbu- 
lent damping respectively) and plasma parameters - most 
notably the CR number density. The streaming speed is 
thus a function of both position and time, and is best self- 
consistently solved in a time-dependent calculation, as we 
will soon tackle]^ Before we forge ahead and use these ex- 
pressions, there are several potential complications worth 
discussing. 

Our streaming speeds for the ICM are characteristi- 
cally of order the Alfven speed, although this can vary spa- 
tially and temporally as plasma parameters vary, particu- 
larly the CR number density. 



Enl31in et al. (20111 argue 



against the Alfven speed as a characteristic CR propaga- 
tion speed in a high /3 plasma, arguing that in the limit 
where the background magnetic field B — >■ 0, this would im- 
ply that v-D VA ^ 0, rather than vu — >■ c, as might be 
expected if there is no magnetic field to couple the CRs to 
the plasma. Instead, they advocate the sound speed cs as 
a characteristic streaming speed. We have several remarks. 
The streaming speeds we have calculated via quasi-linear 
theory assumes {5B/B) <^ 1, and we have checked that this 
condition is self-consistently fulfilled in the ICM (typically, 
{5B/B) ~ 10"'*), an amplitude similar to that inferred for 
the coronal gas in our Galaxy. The hypothetical limit B — >■ 
(which is not realized in the ICM) clearly violates this as- 
sumption, and requires a fully non-linear calculation. There, 



A||(Ax)/Ax 2> 1, whereas we have (A||/Ax)min < 1- Typical ed- 
dies in the MHD cascade vary mostly in the transverse direction, 
k± S> fcy , whereas we seek waves injected by CRs with the least 
possible transverse variation, fex fcy ■ 

^ We have also assume n(> 7) to be a fixed power-law, whereas 
it steepens with time due to energy dependent streaming. 
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we might expect that instabihties generated by a current of 
streaming CRs (e.g., ^Bell (1978) ) would nonetheless gener- 
ate a B-field which will confine the CRs. Nothing in our 
calculations singles out the sound speed as a reference ve- 
locity. 

The resonance condition, equation ([T]), shows that CRs 
of larger pitch angle (^ — >■ 0) interact with waves of progres- 
sively shorter wavelength. However, growth rates Fcr oc /i 
(e.g., Kulsrud (20051), while non- linear Landau damping 
Fnl oc l//i (equation ([3|), so there is relatively little en- 
ergy in such short wavelength waves as — > 0. On the face 
of it, this would imply that it is impossible for particles 
to scatter across the 9 — 90° point via resonant scatter- 
ing to reverse direction, the well-known '90° problem' (in 
fact, the affected region is small; quasi-linear interactions 
can effectively scatter CRs down to /ic ~ 10~^. The gap is 
a little larger, ~ Ui/c ~ 3 x 10~^, in a high-/? plasma when 



ion-cyclotron damping is effective (Holman, lonson & Scott 
[1979] )). The fact that CRs appear to be efficiently confined 
and isotropized in our Galaxy implies that Nature has found 
a way around it. The leading explanation appears to be mir- 
ror interactions from MHD waves created by the 6 ~ 0° CRs, 
which are able to trap the particles and turn them around 
I Felice fc Kulsrud|2001| |. These mirror interactions can also 
be thought of as resonance broadening ( [Achterberg |1981 ' 
Yan fc Lazarian||2008 1 of the long wavelength waves. |Felice 



& Kulsrud (20011 conduct a detailed boundary layer calcu- 



lation of the mirror interaction and find that it introduces 
a minor logarithmic correction (which we have ignored) to 
the standard calculation. We note that if there were indeed 
a 90° problem in the ICM, the resulting light-speed stream- 
ing speeds would imply flat CRp and CRe profiles, which 
is inconsistent at least with observations of radio mini-halos 
such as Perseus. It would also shut ofl radio halos extremely 
rapidly, regardless of how the relativistic electrons are pro- 
duced. 

The wave damping rate is the sum of all damping pro- 
cesses, and thus in principle one should always consider the 
contribution from both turbulent and non-linear Landau 
damping. In practice, we consider limiting regimes where 
one process dominates. Their ratio is: 



Ftur 



3/2„l/4 rl/2 



■ rl/2 ^1/4 1/2 

^MHD,100-'4kcV"CR,-10 



n/2-2 
7l00 °C 



1/2 



(16) 

Turbulent damping thus always dominates at late times as 
the CR profile falls (ncR — ^ 0) and flattens (Lz — >■ oo). 

We have only considered CR self-confinement, and ig- 
nored other possible mechanisms for scattering CRs. As 
we have previously discussed, the anisotropic nature na- 
ture of Alfvenic MHD turbulence (which is mostly trans- 
verse on small scales comparable to the gyro-radius, in con- 
trast to the parallel modes required to scatter CRs) make 
them inefficient scatterers of CRs ( Chandran||2000 Yan & 
|Laza rian"2004l). While the distribution of slow magnetosonic 
waves follows that of Alfven waves (Lithwick & Goldreich 



20011, fast magnetosonic modes can potentially have an in- 



dependent non-linear cascade which is isotropic and can effi- 
ciently scatter CRs ( Schlickeiser|2002 Brunetti & Lazarian 
2007). For now, we eschew this possibility, in favor of the 



well-established self-confinement picture, which is the gen- 
erally accepted theory in our Galaxy. One failing of the self- 



confinement picture in our Galaxy is that both non-linear 
Landau damping and turbulent damping appear to damp 
the waves too efficiently at high energies; the increase of 
streaming speeds with energy appear inconsistent with the 



low observed CR anistropy for i? > 100 GeV (Farmer & Gol 
|drcich "2004j). Chandran ('2000') has proposed that magnetic 



mirror interactions in dense molecular clouds could provide 
this further confinement, though the possibility remains that 
some aspects of the physics are still not well understood. A 
conservative reading of these possible complications would 
take our derived streaming speeds and diffusion coefficients 
as upper bounds; they could potentially be lower if scatter- 
ing is more efficient. 



2.2 Cosmic-Ray Transport 

2.2.1 Cosmic-Ray Transport Equation 

The cosmic ray transport equation in the limit of large wave- 
particle scattering is ( |Skilling|1971| ): 



dt 



(u + va) ■ V/p = V • (Kpnn ■ V/p) 



(17) 



Here, /p(x,p, t) is the cosmic ray distribution function 
(isotropic in momentum space), u is the gas velocity, va 
is the Alfven velocity, n is a unit vector pointing along 
the magnetic field, and Q is a cosmic ray source function. 
Throughout this paper, we shall always use the 3D distri- 
bution function /p, which does not include the differential 
volume factor 4Tvp^ . All momenta p, unless otherwise speci- 
fied, will always be in units of mc throughout this paper. The 
actual momentum will be written pi — pirriiC, the subscript 
denoting the particle type. Any distribution functions writ- 
ten as functions of particle energy rather than momentum 
will be related by 



dui = Atvp'^ fi{pi)dpi = fi{Ei)dEi 



(18) 



E, = 



2 2 

- pfniiC 



dE, = 



PiUiiC^dpi 



Equation (17 1 is derived from the collisionless Vlasov equa- 



tion, which expresses conservation of phase space density: 



dl 
dt 



+ V • (/v) + • / 



dt 







(19) 



but evaluated in the frame of the Alfven waves (which has 
velocity u + va, the sum of the local gas and Alfven ve- 
locities). The distribution function is then expanded in in- 
verse powers of the CR-wave collision frequency I'^f ~ 
/o + /i + /2 + •-, where fr = 0{i^~''). Equation ml is 



obtained after averaging over pitch angle (justified in the 
limit of frequent scattering), and is accurate to second or- 
der, 0{i'~^). The term with Hp expresses diffusion relative 
to the wave frame, and is discussed in detail below. For 
the details of this expansion we refer the reader to |Skillingl 
(19711. This equation implicitly assumes that fo ^ /i, i.e. 



to leading order strong wave-particle scattering renders the 
distribution function isotropic in the wave frame. As we have 
seen, for most plasma parameters (t^D — va)/c <^ 1, so this 
assumption is justified. 
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The physical interpretation of equation (17 1 is easy to 
understand. The left-hand side of this equation is a total 
time derivative, including an advection term in the frame of 
the waves. The first two terms on the right-hand side repre- 
sent diffusion along magnetic field lines relative to the wave 
frame and adiabatic losses/gains respectively. As long as we 
have a functional form for Kp and Q, this equation com- 
pletely describes the evolution of the cosmic ray population. 

Note that in the frame of the wave, and considering the 
isotropic part of the distribution function /o (so that there 
is no diffusion relative to the wave frame, Kp — 0), we have: 

Dfo 



(20) 



with 



i.e. the C Rs evolve adia baticallv in the wave frame, 
p oc (Skilling 19711. This makes physical sense: there 

are no electric fields in the frame of the wave, and hence 
the particles conserve energy; they can only scatter in pitch 
angle. However, the CRs do not evolve adiabatically in the 
frame of the gas, where there are electric fields associated 
with the hydromagnetic waves. Thus, there is an irreversible 
energy transfer from the CRs to the gas, with volumetric 
heating rate (e.g., Kulsrud (2005)): 

(21) 
This 



Twave — — VA ' V-Pc, 

which we shall refer to as the "wave heating rate" 
may be thought of as the rate at which work is done on 
the gas by CR pressure forces, va • F. Importantly, this 
heating rate is not Fwavc = — vd- VPc, as has sometimes been 



adopted elsewhere in the literature (e.g., Uhlig et al. ( 2012 1) 



The latter expression gives rise to unphysically large heating 
rates when vo 3> wa. Super- Alfvenic streaming arises due to 
a reduction in coupling between CRs and gas; it is unphysical 
that this would give rise to greater heating. Physically, all 
momentum and energy transfer between the CRs and gas is 
mediated by hydromagnetic waves; the rate at which work 
is done by any transmitted forces is therefore set by the 
velocity of the waves va. 

To next order in , slippage with respect to the wave 
frame is expressed by the diffusion coefficient Kp-. 



k(7) 



2^(^,7) 



(22) 



where the wave-particle collision frequency i^(/i, 7) is (Kul 
|srud fc Pearce|[Tt69l ): 

l3 



(m,7) 



(23) 



and the average is taken over pitch angle (Skilling 19711. 
This expression is obtained from equation (|4|| as shown 
by Kulsrud & Pearce (19691, and we assume relativistic 



CRs such that i; ~ c. From equation (221, the more fre- 



quently CRs interact with Alfven waves, the more slowly 
they diffuse relative to the waves — as one might expect, 
since scattering isotropizes the CR in the wave frame. Equa- 
tion ( |23[ ) can be understood from the fact that a single CR- 
wave encounter in one gyro-period r leads to a change in 
pitch angle M ^ {SB/B) (|Kulsrud|[2005|; thus N ~ t/r 



encounters leads to a net random walk in pitch angle of 
{Aef ~ N{5B/Bf ~ t/T{5B/Bf, or a pitch angle diffu- 
sion rate of De ~ {AOf /t ~ no{SB/B)'f\ Equation ([22 1 



The mean free path of a CR is roughly the distance over which 



can be evaluated by equating wave growth and damping 
rates to obtain the amplitude of the waves, {5B/Bf, as for 
instance in equation It can also be intuitively written 
in terms of streaming speeds. From equation (171, we can 



write the net streaming speed (i.e. the frame in which the 
mean CR flux vanishes), as ( [Blandford fc Eichler|[T987l ) : 



V-D 



Mp) [ 



V/p 



(24) 



where n is a unit vector pointing along the magnetic fleld, 
down the CR gradient. The first term effectively corrects 
for the Compton-Getting effect, i.e. the differential Doppler 
shifts of particle energies in transforming from the wave to 
the inertial frame (depending on whether particles are mov- 
ing parallel or anti-parallel to the wave, when we calculate 
the particle flux in the inertial frame, we must compare par- 
ticles of slightly different energy in the wave frame). If we 
solve this for the diffusion coefficient, we obtain: 



•V/p 



-vu — 



1 dloi 
3 91ogp 



/p 



L^[vo--va] (25) 



Here we have set 91og/p/91ogp « —4.5; as before, the en- 
ergy dependent CR scale length is I/z(7) ~ j/p/(n ■ V/p)|. 



If we insert this into the diffusion term in equation (171, we 
obtain: 

D{r) = V • (Kpnn • V/p) ^ V • (/pn(t;D - va))- (26) 

where we evaluate the drift speed relative to the wave frame, 
{vu — Va), from equations (10 1 and (15 1. Note that the 



gradient of the distribution function V/p (or equivalently, 
the scale height Lz) does not appear in the diffusion term. 
It only appears if (wd — ua), has a functional dependence 
on Lz. This is true for non-linear Landau damping, where 
(I'D — '"a) oc L7^^^ (so that the diffusion term oc (V/p)^^^ 
rather than (V/p)), but false for turbulent damping, where 
the diffusion term is therefore independent of the magnitude 
of V/. 

The latter unusual behavior was first noted by [Skillingl 



( 1971 1 for the case of ambipolar damping, which shares sim- 



ilarities with turbulent damping in this regard (although 
he dismissed it as unimportant, since the effects of diffu- 
sion were small for the applications he considered). From 



equations (22 1 & (231, and equating wave growth and with 



a generic damping rate Fd, the diffusion term can be ex- 
pressed more transparently as ( Skilling|1971 1: 

52, 



D(r) 



/ Fd-B^ii n • V/p 



yA-K^nipQoVA 
1 „ 



47r3p7/2gl/2j„^. 



1/2 



L-V/p 

B^n. 

(-1/2 In 



V/p 
V/p I 



(27) 



(28) 



where we specialize to the case of turbulent damping in the 
second equality, and substitute Fturb ~ '^A/('-LiMHD)"^'''^. 
Note that, other than the sign of n ■ V/p, the term within 
the divergence is independent of fp. This has important con- 
sequences for us, in that diffusion does not slow down with 



the pitch angle diffuses by order unity (so that the CR reverses 
direction), A ~ cDg^ --^ 3 X 10^'^(5B/B)~^ cm, where the pitch 
angle diffusion coefficient Dg ^ {SB / B)~'^Clo. Thus, even small 
fields of (SB/ B) ~ 10~^ would lead to mean free paths of A ~ 1 
pc, implying that the diffusive approximation is excellent. 
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time as V/p decreases. Instead, it is independent of /p and 
depends only on plasma properties. If these plasma proper- 
ties are roughly constant over the streaming timescale, then 
fp{r,p,t) « D{r,p) « is roughly constant and tstroam oc 
/p//p oc /p, with decreases with time as fp falls. This accel- 
eration is key in our more detailed calculations which show 
that large changes in radio halo luminosity are possible de- 
spite apparently long initial diffusion times. It is important 
to stress, however, that while the diffusion time with turbu- 
lent damping is not sensitive to the magnitude of V/p, it is 
still sensitive to the sign of V/p. The sign of n • V/p reflects 
the fact that CRs can only stream along B-fields, down their 
gradienlQ diffusion has no further effect if V/p = 0. Failure 
to carefully treat this can result in spurious numerical insta- 
bilities ( [Sharma, Colella fc Martin||2010 1, which we discuss 
in SJSH 



2.3 Collisional Losses 

Cosmic ray protons can also lose energy from direct collisions 
with gas particles, either through Coulomb interactions, or 
hadronic interactions (pion production). While these are 
generally subdominant to losses from wave-particle inter- 
actions, we include them for completeness. This transfer of 
energy from CRs in turn heats the gas. 

The energy loss rate of a CR of speed /3 = v/c and 
kinetic energy E due to Coulomb collisions in ionized gas is: 
pvfannheim fc Schlickeiser| ( [1994] )) 



-4.96 X 10-i^erg s'^ ( 

V cm / 



(29) 



Here a;^ = 0.0286[r/(2 x lO'' K)]^/^ with T and the 
gas electron temperature and number density. The energy 



loss rate of a CR due to hadronic collisions is ( Mannheim & 



Schlickeiser (1994 1) 



cLE 
d^ 



'.0.5nNapp/3cEe{E - Etbr) 



(30) 



where the pp cross section for hadronic interactions is o-pp 
and the target nucleon density is riN = We/(1 — 0.5Y), Y 
being the helium mass fraction. The above assumes an in- 
elasticity of A' = 1/2 for the collision. The Heaviside step 
function enforces the condition that only cosmic rays with 
kinetic energy above i5thr = 282 MeV undergo pion produc- 
tion. All of the energy loss in Coulomb collisions goes toward 
heating the gas, whereas only ~ 1/6 of the inelastic energy 
in hadronic collisions goes toward secondary electrons which 
heat the gas, the rest escaping as gamma rays and neutrinos. 

These loss terms are represented in the CR Vlasov equa- 
tion as: 



dt 



C,h 



= -|(pc,h/p) 



(31) 



PC,h = 



^ CRs can only stream up a gradient if the sign of energy transfer 
is reversed - i.e., the gas gives energy to the CRs, rather than 
vice-versa, as in Fermi acceleration. In this case, the picture of 
self-confinement is clearly not applicable. 



where E = (x/lT+l? — l)mpC^ and the momentum p is in 



units of rUpC. 



2.4 Turbulent Diffusion 



As we have seen, turbulent gas motions can damp MHD 
waves and enhance CR streaming. However, they can also 
directly transport CRs advectively. A proper treatment of 
the interplay between these effects requires 3D MHD sim- 
ulations. Here, we will simply treat turbulent motions as a 
diffusive term in the CR transport equation. If PcR/^gas 
is small and CRs have negligible effect on the dynamics, 
they simply act as a passive tracer species. Analogously to 



the mixing of metals by turbulent diffusion ( Rebusco et al 



20061, we can write: 



dncF 



-V ■ 



tturb"-EV 



where 



vtLt vaL^ 



(32) 



(33) 



i.e., if turbulent mixing is vigorous, the CRs will have uni- 
form relative abundance, ncn oc n^. This has some support 
from simulations where CR dynamics are taken into account 
( [Sharma et al.|[2009[ ). There, turbulent convection results 
in constant CR entropy -P/^^c^^ (where 7cr = 4/3) and 
Pen/ Pg, =const. This implies ncR oc P^^'"^^- . Since stratified 
gas in a cluster has a polytropic equation of state Pg oc Pg*"' 
where 7pt ~ 1.2 — 1.3 (e.g., Capelo, Coppi fc Natarajan 



(|2012|), this implies ncR oc (p^)Tpt/-TCR ^ p^ "-" "", con- 
sistcnt with our assumptions. Alternatively, [Enfilin et al.| 
(2011) suggest a target profile set by gas entropy, rather 
than CR entropy: ncR oc Pg^"*" , where 7g = 5/3. In this 



case, all occurrences of ne{r) in equation (321 will be re- 



placed by ri{r) 



and ncR oc pl^^^''^ 



Given the many uncertainties in the model, this difference 
in scalings is of secondary importance. 

We also need to take adiabatic heating and cooling into 
account. The normalization of the distributio n function f 



i/3 



(e.! 



varies with adiabatic changes as C oc n^j^ 
( 2007 1), where a = 4 — 5 is the spectral slope of the distribu- 



Enfilin et al. 



tion function. Thus, the overall efi'ect of turbulent diffusion 
on the distribution function is: 



dl 
dt 



= -V 



Kturb5"''^V ^ 



/ 



(34) 



and S{r) = P^ 



1/7CR 



,{r), or 5{r) = ri{r) = P< 



l/7g 



These equations show that turbulent diffusion, acting 
alone, will lead to a centrally peaked CR profile similar to 
the gas profile. On the other hand, turbulence also damps 
MHD waves, leading to enhanced outward streaming, which 
flattens the CR profile. Which effect dominates? While we 
explore this in detail in our numerical calculations, it is use- 
ful to first get an order of magnitude estimate. From equa- 



tions (25 1 and (33 1, we obtain: 



'^'Stream 



VA 



^MHD 



r3/2 



(35) 



We expect Lz/Lmhd <; 1, and {vo/vA— 1) ^ 1 (from equa- 
tion (15l) for turbulent damping; moreover, these factors 



© 0000 RAS, MNRAS 000, 000-000 



Cosmic Ray Streaming in Clusters of Calaxies 9 



increase during the streaming processes as Lz rises and ncR 
fall. Thus, 

/^stream <; Kturb in our fiducial model. Moreover, 
if the strength of turbulence increases such that Lmhd falls, 
^stream /«turb rises. Stronger turbulence has a larger effect on 
damping of MHD waves than on inward advection of CRs, 
and the CRs stream outward faster. Thus, in this framework, 
turbulent diffusion can never establish a centrally peaked 
profile, regardless of its strength. In practice, coherent bulk 
motions (triggered by mergers, or perhaps by gas 'sloshing') 
can potentially bring CRs to the cluster center, and/or pro- 
duce a magnetic topology which is unfavorable for outward 
streaming. However, modeling such stochastic events is be- 
yond the scope of this paper. 

A few comments about our choice of fiducial parameters 
for turbulence is in order. It is customary to define [Lt,vt), 
where Lt is the outer scale, and Vt is the velocity at this scale. 
Instead, we work with (Lmhd, «a), where Lmhd is defined to 
be the scale at which the turbulent velocity is va - In general, 
Vt ~ va{Li/ LmhoY^^ , and more vigorous turbulence can be 
characterized by smaller values of Lmhd. However, if there is 
equipartition between = B'^/8n and Ut = l/2pVt, then 
Vt ~ va and thus Lt ~ Lmhd- Thus, (Lmhd, va) are sensible 
fiducial parameters. Secondly, we have assumed that Lmhd 
(or equivalently, Vt at a fixed scale) is independent of ra- 
dius. Is this consistent with cosmological simulations, which 
show that turbulent pressure support becomes increasingly 
important with radius? A fit to low-redshift clusters gives 



spherically symmetric version of ZEUS3D, previously used 



( Shaw et al. (20101, see also Battagha et al. (2012l) 



Lturb 
Lthcrm 



— ao 



Rs 



(36) 



where ao ^ 0.18 ± 0.06 and Unt = 0.8 ± 0.25. This implies 
Vt oc 7-0 '*~"r/2^ where T oc r~°"^ , and qt is ge nerally small 
(e.g., QfT ~ x^/(l + 1.5a;) (Loken et al. 20021, where x = 
r/rvir, so qt ~ 0.2 at r = 0.5rvir). On the other hand, 
given our assumption that B oc p"^ (see S3.4I, we have 
VA oc r'"'^""^'"", where p oc r'"'" , and Qp ^ 2 — 3 over 
most of the cluster. The radial scalings for vt and va are 
thus roughly consistent: for instance, qb ~ 0.3 (as assumed 
for Perseus & Coma) gives va oc j-O ^-O-S Finally, we note 
that for our assumed levels of turbulence, heat dissipation 
is relatively unimportant. The heating time is: 



[Ahcrm / L'thc: 

f'hcat ~ ^ f'turb 



L'turb 



5Gyr 



.A.sLk 



''A, 100 



(37) 

where ft,5 = [(fthcrm/t^turb)/5], Lmhd, 100 = (Lmhd/100 
kpc), and i^a.ioo = (t'^/lOOkms""'^). 



3 METHOD 

Our main task is to solve the CR transport equation, in the 
form: 



dt 



+ (u + va) ■ V/p = V • (Kpnn • V/p) 



+ ^pfv.(u + VA) + Q-| 



(PC.h/p) 



(38) 



where the last two terms are as in equations (31 1 and (341 



to solve the CR equations in the fluid approximation (Guo 
irOhl[2008t 



Our goal in this paper is to determine if CR stream- 
ing is a plausible means of turning off radio halos in the 
hadronic scenario. We therefore run numerical simulations 
where the cluster is assumed to be in strict hydrostatic and 
thermal equilibrium, and only solve the CR transport equa- 
tion (ignoring the fluid equations for the gas, equations , 
by setting all time derivatives to zero) to examine the ef- 
fects of CR streaming. In the absence of a cooling flow, the 
only time-dependent terms in the fluid equations for the gas 



(equations (44 1, (45 1, and (46l) relate to the CRs, and have 
negligible effect. We initialize the CR profile so as to repro- 
duce the observed radio surface brightness profiles in the 
classical hadronic model, and follow the time evolution of 
radio emission as the CRs stream out. In this methods sec- 
tion, we discuss numerical regularization of CR streaming 
(j ]3.1| l, a test comparison of our CR transport solver in the 
fluid approximation {[3.2 \, calculating radio and gamma- 



ray emission (^3.31, and our initial conditions (S3.4I for a 
prototypical radio mini-halo (Perseus) and a prototypical 
giant radio halo (Coma). Results are then presented in the 
following section, ij4] 



3.1 CR streaming: Numerical Stability 

Cosmic rays can only stream down their gradient, in a di- 
rection: 



s = -sgn(B ■ V/p) 



B 

IbT 



(39) 



In our ID simulations, s — — f sgn(d/p/dr). However, if this 



is enforced in equation ( 28 1 , it leads to numerical instabili- 



ties and unphysical oscillations in the distribution function. 



The origin of this difficulty is easy to understand ( Sharma, 



Colella & Martin 2010 1; it essentially arises at local extrema. 



If the simulation at any time produces a local density max- 
imum, diffusion out of the local maximum will cause the 
density to drop significantly there. If the time step is not 
properly restricted, this decrease will overshoot, causing the 
density to drop below neighboring regions, creating a local 
minimum. The opposite problem will then occur, with in- 
wardly diffusing CRs causing the CR density to increase too 
much. The result is an unphysical oscillation that eventually 
spreads out to all space. Because CR streaming results in a 
flat proflle where (d//dr) vanishes everywhere, this problem 
can become acute as time goes on. 



Sharma, Colella & Martin (20101 show that for an ex- 



plicit code (such as ZEUS3D), the restriction on the time- 
step such that new local extrema are not created is: 



At \ f"\Axyf\v\ 



(40) 



respectively. To do so, we have written a new module in a ID 



which is much more onerous than the standard Courant con- 
dition (we have explicitly verifled that simulations which sat- 
isfy the Courant condition suffer from spurious oscillations). 
They suggest regularizing the CR transport equation by re- 
placing the discontinuous sgn(/p) with the smooth function 
tanh(/p/e) for some choice of e. As e tends to zero, the tanh 
function approaches the sign function. This effectively sets 
e as a minimum scale value for fp-. if /p <C e, the simulation 
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behaves as if /p = 0, and suppresses CR streaming. Alterna- 
tively, it can be viewed as introducing a difTusive term at an 
extremum, with diffusion coefficient /p/e, similar to the use 
of explicit viscosity to regularize Euler/Burger's equations. 
In this case, the maximum time step allowed to suppress the 
instability is: 



(41) 



For us, the relevant speed v is the streaming speed, calcu- 



lated via equation (24 1, which we insert into this equation. 



We have found that a scale value of: 
£ = /p/L; L = 3Mpc 



(42) 



is sufficiently small that decreasing it any further does not 
significantly change the results. In Fig [T] we show a con- 
vergence test (showing the radio luminosity of Perseus as a 
function of time when Lmhd = 100 kpc; see Figjsjs) where 
the figure converges to the correct solution as e is decreased; 
a value of e half of our fiducial value (e = 10~^^/p) gives 
identical results. 

In practice, although we use a smoothing scale e, we use 
the time constraint ( 40 1 rather than ( 41 1 and we use v = va 
rather than v = Vs- Additionally, to prevent the time step 
from dropping to zero we impose a minimum time step 



At ^ 1 X 10"^A2;^e/2/pl?;| 



(43) 



where the numerical factor out front is arbitrary. We do all 
of this to regulate the runtime of the simulation - the less 
stringent time steps will be less accurate but will run quicker. 



and we can adjust the minimum time step ( 43 1 to the desired 



balance of speed and accuracy. As a result, local minima 
and maxima do develop at some points of our simulations, 
however this will always happen as the CR profile flattens; 
as long as the local extrema do not grow unstably the results 
should be robust. 

In addition, we enforce the constraint that A/p ^ 
0.05/p in a single time-step (and similarly for the gas en- 
ergy and density) . Note that while ( 40 1 is only applied 



at local extrema, this condition is held everywhere. Even- 
tually, as /p falls, this shrinks the time-step to zero. To 
avoid this, we define a momentum dependent minimum 
/p,min(p) = 10"^/p(rmax,p, to), where Tmax is the outer 
boundary of the simulation, and to is the initial time. The 
distribution function /p is then never allowed to drop below 
this value. Also, once /p falls below 50/p,min(p) anywhere, 
all time step restrictions there are ignored, including ( 40 1 . 



3.2 Test Case: AGN Feedback 

To test our solver for the CR transport equation, it is useful 
to compare against previous results where CRs are treated 
in the fiuid approximation. Specifically, we compare against 
the results of |Guo &: Oh| ( [2"008[ ), which simulates the effects 
of CRs injected by a central AGN on the thermal state of a 
cool core cluster. It was found that a combination of electron 
thermal conduction (at some fraction / of the Spitzer value) 
and CR mediated wave heating was sufficient to stem a cool- 
ing fiow. The following governing equations for the two-fluid 
ICM (gas and cosmic rays) were used: 



|+V.(pu)=0 



(44) 



£=1x10-24 f 
£ = 3x10-25 
£=1x10-25 

5xio-2Ei r 




3 4 

t (108 yj.) 



Figure 1. Convergence test for the smoothing scale parameter e. 
We plot the 1.4 GHz radio luminosity of Perseus, which declines 
with time due to cosmic-ray streaming (here we assume -Lmhd = 
400kpc; see j]4.1| for details). As e is decreased, our calculations 



converge. Our fiducial value is e = 10 /p. 



dt 



+ V • (Sw) = -VPg - VPc - pV$ 



dt 



+ V ■ (Bgu) = -PgV • u - V ■ F 



dt 



nlK{T) + f],n,E, - va • VPc 



(7c - l)(u -f va)-VSc - V ■ Fc ( 



Fc = 7c-Bc(u + va) - n«:c(n • V 



(45) 

(46) 

(47) 
(48) 



where p is the gas density, Pg is the gas pressure, iJg is the 
gas energy density, S = pu is the gas momentum vector, Ec 
is the cosmic ray energy, Pc = (7c — 1)-Ec is the cosmic-ray 
pressure, F is the electron conduction heat flux, and $ is 
the gravitational potential. The term rjcUcEc, where r/c = 
2.63 x 10~^® cm'^ s~^, takes Coulomb and hadronic heating 
of the gas by cosmic rays into account. The initial conditions, 
gravitational potential $(r), and cooling function A(T) are 
as spelled out in |Guo &: Oh| ( |2008[ ); please refer to the paper 
for details. The source function Q represents the injection 
of CRs by an AGN, triggered by gas cooling: 



Qc 



(49) 



Here, e = 3 x 10 is an efficiency parameter, u — 0.3, and 
ro = 20 kpc is a scale distance. 



The code of Guo & Oh ( 2008 1 uses the CR energy den 



sity Ec as the fundamental dynamic variable for CRs. It is: 



E, 



P Tp{pp)fp{pp)dp 



where 



Tp{p) 



(50) 



(51) 
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is the kinetic energy of a CR proton of momentum pp. By 
using Ec as the main CR dynamic variable, all momentum 
dependence has been integrated out. By contrast, we wish 
to retain momentum dependence, and instead use f{r,p,t) 
as our fundamental variable. We therefore continue to solve 
equations (441 — (46 1, but replace equations (47 1 & (481 



with the equation for the distribution function, equation 
(381, and solve for Ec as required in equations (45 1, (461 



via equations (50 1 and (51 1. In calculating the momentum- 
dependent source function Q for use in equation ( 38 1 , we 
suppose that the injected spectra has the form: 



ME) = 



Ac.eiE - El] 



{E/E,)^ + {E/E,)^^- 



2)/2 



(52) 



which assumes a steady state spectrum at low (high) ener- 
gies due to Coulomb (hadronic) losses, and smoothly con- 
nects these regimes (see ( |Guo fc Oh|2008| for details). Here, 
-E* = 706 MeV is a cross-over energy separating the low- 
and high-energy regimes, while the assumed spectral index 
is a = 2.5, and cutoff energy is E\ = 10 MeV. We create 
a source function with the same momentum dependence as 



equation (52 I, and then normalize it to the total CR injec- 
tion rate given by equation ( |49[ ): 



> 



New code 
Old code 




0.2 



0.4 0.6 

t/t„ 



0.8 



Figure 2. Temperature versus simulation time at some select 
radii for the new code, where we solve the CR transport equation 
( |38| l. The results of the old code, which treats CRs in the fluid 
limit (equations ( |47| l, ( |48[ l), are displayed in the dotted lines. 



=47r / p TpQdp. 



(53) 



To maintain consistency with Guo & Oh (20081, we use 



the same momentum-independent diffusion coefficient used 
there. 

The simulation grid has two ghost zones at each end in 
the radial direction, and one ghost zone at each end in the 
momentum direction. The density and temperature of the 
ICM are linearly extrapolated into the spatial ghost zones. 
For the CR spectrum, constant boundary conditions are en- 
forced in the radial direction, i.e. /p(p) at the spatial ghost 
zones are set equal to /p(p) in the adjacent active zones. In 
the momentum direction we required that dlog/p/dlogp be 
constant across the boundary. As for time step constraints, 
for this test case we do not allow pg or Eg to change by more 
than 25% in any time step. We also enforce the Courant con- 
dition for all cells. At < Ax^/2Kp. 

We find that our full model reproduces the results of 
|Guo &: Oh| ( |2008[ ) extremely well. A example is shown in 
Fig. [2j where we show the temperature as a function of 
time for several select radii. The cluster is initialized to be 
isothermal; after an initial transient, it is thermally stabi- 
lized against a cooling catastrophe by a combination of CR 
heating and electron thermal conduction. The dotted lines 
show the results from the code of |Guo fc Qh| p008), which 
integrates the fluid equations, while the solid lines indicate 
the results of the new code, which computes the distribution 
function. Indeed, even when we include the full momentum 
dependence of the diffusion coefficient, the results barely 
change (for this particular example, we have assumed non- 
linear Landau damping) . This is because most of the energy 
density of CRs is dominated by low energy CRs (~ 1 GeV) 
for which the diffusion time is negligibly long. As we shall 
soon see, diffusion cannot be neglected for the high-energy 
CRs which are responsible for observed radio emission. 



3.3 Computing Emissivities 

Our fundamental simulation variable is the CRp distribu- 
tion function, fp{r,p,t). Here, we describe how radio and 
gamma-ray emission can be inferred from /p(p) (hereafter, 
we suppress r, t) in the hadronic model, given an assumed 
gas density and magnetic field. 

How is radio emission produced? CR protons undergo 
hadronic interactions to produce pions, which in turn de- 
cay to produce relativistic electrons (CRp + nucleon — > 



27). 



The high energy electrons which produce observable syn- 
chrotron emission have short cooling lifetimes and we there- 
fore assume a steady-state between injection and cooling. 
A CRp distribution /p(p) gives rise to a pion source func- 
tion due to the hadronic pp interaction (Pfrommer, EnfJlin 
fc Springel|2008| ): 



dPp/p(Pp)cnN^(Pp)o-p 



X 5 ( P. - ^Pp ) e(pp - 0.78) 



(54) 



where crjp = 32(0.96 + e^-'^''^-'^"^) mbarnR The delta func- 
tion enforces the mean pion momentum (pT) ~ Pp/4 and the 
Heaviside step function 9 incorporates the threshold proton 
momentum for the pion production to occur. Approximating 



We follow [Pfrommer fc Enfilin| |2004[ | in absorbing the weak 
energy dependencies of the pion multiplicity and the inelastic 
cross-section in this semi-analytical parametrization of the cross- 
section, where Op is the average CR spectral index. 
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the pion multiplicity ^ as 2, this gives: 
, , 16 TT . / 4m^ 



3 mn 



4m^ 



(55) 



X /p ( ^^P^ 1 e ( ^^p^ - 0.78 
V mp y \ 'Tip 

Under this approximation, the neutral pion source function 
s^o is the same. 

The charged pion population will undergo pion decay, 
producing electrons (and other particles). This electron pro- 
duction is described with the electron source function: 

Ap-n 



Sc{pc) = S^±(p^(Po)) 



dpe 



, , 64 mc ^ 

So(Po) criNCr 47r 

3 rrio 



16mc 



X /p ( Pvr 1 e ( ^:^Po 



0.78 



(56) 



(57) 



In the second equation we have used — 4pe. If we as- 
sume an equilibrium between this source and any losses, i.e. 
a steady state solution, the electron spectrum is then deter- 
mined from 



.fe(Pe) = 

where the losses pc are 

Pe(pc) 



Ap'cScip'o] 



4 (TTCPc 



3 nicC^ 



(eb + £c 



(58) 



(59) 



from synchrotron radiation and inverse Compton (IC) scat- 
tering. Here, (tt is the Thompson scattering cross section, 
and eb , Ccmb are the energy density of the B-field and cosmic 
microwave background. 

From the electron distribution function we can deter- 
mine the resulting synchrotron emissivity ( [Rybicki fc Light^ 
man||1979|): 



>(r) = 0.333 



2n rricC^ 



d7e.fo(r,7o)F 



(60) 



In the above, Vc = 3eB7c/47rmcC and the function F 
is an integral of a modified Bessel function, F{x) — 
X K^^3{x')dx' . The numerical factor in front comes from 
averaging the CRe population over pitch angle, assuming 
isotropy. The observed surface brightness is: 



5.(ri) 



(1 + ^) 

where 1^' — v(l + z). The luminosity is 



3^' (r) 



r dr 



(61) 



(62) 



We also make predictions for gamma-ray emission. We 
only consider gamma-ray emission from neutral pion decay 
7r° — ^ 27 and ignore the subdominant contribution from in- 
verse Compton scattering. The analysis is much the same as 
above. Following Mannheim & Schlickeiser ( 1994 1, we derive 



a photon source function from the pions: 

dEnS^oiE^:) 



(63) 



where the neutral pion source function s^o(E-n) is assumed 



to be the same as for charged pions, equation ( 55 1. From this 



source function, we determine a number production rate per 
unit volume A^.: 

/•oo 

\-,{>E^)=\ dE\s^,(E\) (64) 
and the flux detected at Earth above an energy E'^ \ 



- 4^ 



d^rA.^(> E^) 



(65) 



Given the strong momentum dependence of CR stream- 
ing, it is worth clarifying which range of CRp momenta are 
most observationally relevant. For radio emission, the char- 
acteristic synchrotron frequency is ~ 37'^i/c, where Vc is the 
non-relativistic synchrotron frequency. For a given observa- 
tional frequency Va, the greatest contribution comes from 
electrons with: 



Pcmit ~ 7cmit « 4 X 10^ ( 1 GHz ) 



1/2 



B 
3^iG 



-1/2 



(66) 



Thus, ~ 10 GeV electrons are responsible for ~GHz emis- 
sion in ^G fields. Typically, pe ~ (l/16)pp, where pe is the 
momentum of a secondary CRe produced hadronically. The 
reduction in energy by a factor of ~ 16 comes from the fact 
that the limiting inelasticity is ~ 1/2 ( Mannheim & Schlick- 



eiser ] p94l ), the pion multiplicity is a factor of ~ 2 due to 



2 pion jets leaving the interaction site (Nachtmann 19901, 
and {E) = (l/4){£;^±) in the reaction tt^ -s- -I- 'iv. Thus 
for ~ GHz emission, ~ 100 GeV protons are most relevant, 
while for LOFAR observations at ~ 100 MHz, ~ 10 GeV pro- 
tons are most relevant. For 7-ray emission, E^, ~ ( 1/8) i5p 
(all the factors are as before, except E-^ = l/2£'.^o). Thus, 
Fermi, which is most sensitive in the E.^ « 0.1 — 3 GeV range 
(rather than 0.1 — 300 GeV, due to the pion bump), probes 
7?p ~ 1 — 30 GeV, while imaging air Cerenkov telescopes 
such as MAGIC, HESS and VERITAS are most sensitive in 
the _E.y ~ 0.3 - 1 TeV range (rather than 0.3-10 TeV, due to 
the steep CRp spectrum), probes i7p ~ 3 — 10 TeV. 



3.4 Initial and Boundary Conditions 

We choose to simulate a prototypical radio mini-halo, 
Perseus, and a prototypical giant radio halo, Coma. We 
choose initial conditions which reproduce current observa- 
tions of their radio surface brightness, and then watch how 
this evolves under the influence of streaming. For the Perseus 
cluster, we adopt empirical fits to the cluster temperature 



and electron density profiles (Pinzke & Pfrommer (2010l) 



based on observed X-ray emission (Churazov et al. (20031) 



10-3cm-3 



46 



1 + 



57 kpc 



-F4.79 



1 + 



200 kpc 



(67) 



r = 7 keV 



1 + (r/71 kpc)^ 
2.3 + (r/71 kpc)» 



1-h 



380 kpc 



(68) 
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From these we determine an internal energy distribution (via 
the ideal gas law) and a gravitational potential (via hydro- 
static equilibrium). Similarly, for Coma the fits are (Pinzke 
& Pfrommer (20101, based on Brief, ffenry & Boehringer 



(19921) 



no 
lO-^ci 



= 3.4 



T = 8.25 keV 



1 + 



1 + 



294 kpc 



460 kpc 



(69) 



(70) 



The radio surface brightness profiles at f .4 GHz are fit by a 
P profile: 

S(r) = So [l + (r/rc)'] (71) 
which is reproduced by an initial emissivity of 

Mr)^j,,o[l + ir/r,f]-'^ (72) 



>,o=2^(6/?-l)^Q,3/3 



where B is the beta function. For Perseus, /3 — 0.55, 



30 kpc, and So = 2.3 x 10"^ Jy arcmin (^Pedlar et al 



1990), while for Coma P = 0.78, rc = 450 kpc, and So = 



1.1 X 10"^ Jy arcmin"^ ( [Peiss et al.|l997 l 



We assume that the magnetic field scales with gas den- 



sity: 



B^Bo 



Ticjr) 



(73) 



Such a scaling is motivated by simulations (Dubois 
& Teyssier 2008 1 and Faraday rotation measurements 



(Bonafede et al. 2010 Kuchar & 



Enfilin||201l| ); for in- 
stance, rotation measurements for Coma are well fit by 
qb « 0.3 - 0.7 ( [Bonafede et al.|2010[ ). In the future it would 



be interesting to explore other scalings, if for instance this 



relationship also has temperature dependence (Kunz et al 



2011 1. We find that radio surface brightness profiles can 



be well fit by ob = 0.3 for both clusters, but for Perseus 
Bo = 10 /iG, while for Coma Bo = 5 /iC. We choose a cosmic 
ray distribution function motivated by cosmological hydro- 
dynamic simulations of galaxy clusters where cosmic rays 



are accelerated via diffusive shock acceleration (Pinzke & 
Pfrommer|2"oTol ): 



/p(r,pp)=C(r)^A,p-' 



(74) 



A = (0.767, 0.143, 0.0975) a = (2.55, 2.3, 2.15). (75) 
and the normalization 

Cir) = ~ ^^:"!_;^ + Ccente.. (76) 



1 + 



Note that these simulations do not take into account 
the effects of cosmic ray streaming. The parameters 
Cvir, Cccntor, nrans are thcu choseu such that the model ra- 
dio brightness profile agrees with fits to observations (equa- 
tion (3.4l). For Perseus, if we define C(r) = C{r)nc{r), 
then Cccntcr = 8.3 X 10-^C'vi^ = 7.2 X 10-^ n^ans = 36 
kpc, /3c = 1.0. For Coma, Cccntor = 6 x 10"" cm"^. 



Cvir = 5.2 X 10"" cm-^, rtrans = 55 kpc, /3c = 1.09. The 
initial radio surface brightness profiles derived from these 
parameters are shown in FigjSi and [5^. 

When we solve equation (38|, the simulation grid has 
two ghost zones at each end in the radial direction, and 
two ghost zones at each end in the momentum direction. To 
set values in the ghost zones, we use dlog/p/dlogp =const 
in the momentum direction at both the inner and outer 
boundary, i.e. a power law extrapolation. In the spatial di- 
rection, we use d log fp/d log r =const at the inner boundary. 
The outer boundary requires a little more care, since it can 
fall to extremely low values which result in round-off error; 
also, if the CR gradient goes to zero at the outer simula- 
tion boundary, this artificially suppresses CR streaming. For 
Perseus, we use dlog/p/dlogr =const at the outer bound- 
ary, but subject to the condition that /min ^ ./imax+i 



^+2 



where 



is the in- 



dex of the last active zone, /min(p) = 10 ^/p(r-max, p, io), 
and X = 0.98. For Coma, where the initial profile is al- 
ready extremely flat, we simply adopt fi^^^+i = Xfi^^^, 
/imax+2 = ^/imax+i- To couscrve CRs during the process 
of turbulent advection, we also enforce the CR turbulent 
diffusion flux, defined as Ftmb = K-tmhS'^^^V {fpS^°'''^) (see 
equation ( |34[ |), to be zero at both spatial boundaries. 

For both Perseus and Coma, we use 1.4 GHz data. Note 
that for Coma, which is the most well-studied giant radio 
halo, recent 1.4 GHz and 352 MHz data cannot be reconciled 
by the classical hadronic model with a power-law spectrum 
(Brown & Rudnick 2011), though this conclusion is subject 
to systematic uncertainties in the zero-point of 1.4 GHz data 
( Zandanel, Pfrommer &: Prada|2012 I. We shall also see that 
energy dependence in the streaming speed alters the CR 
distribution function, so that it is no longer a power-law in 
momentum, potentially solving this problem. 



4 RESULTS 

We now show results for a canonical radio mini-halo in a 
cool core cluster (Perseus), and giant radio halo in a non- 
cool core cluster (Coma), starting from the initial conditions 
given in §3.4| We use Perseus to illustrate most of the rele- 
vant physics. Unless otherwise noted, all calculations assume 
^MHD = lOOkpc, where Lmhd is the lengthscale at which 
?;a = '^turb (note from e = v%/Lmhd that smaller values of 
^MHD correspond to more vigorous turbulence). 



4.1 Perseus Cluster 

The initial conditions for Perseus correspond to a CR profile 
where the ratio of CR energy to the thermal gas energy is 
almost constant throughout the cluster, decreasing slightly 
in the outskirts. In Fig [3^ we compare radio emission from 
our initial conditions to surface brightness observations at 
1.4 GHz from [Pedlar et al.| ( [T990l ). Note that the observa- 
tions only span a limited radial range (which produces ~ 1/2 
of the total radio luminosity in our model). The normaliza- 
tion of the profile falls substantially in several hundred Myr, 
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while its shape does not evolve significantlj]^ Figjsjj shows 
the evolution of the 1.4 GHz radio luminosity with time, 
and how it depends on the strength and nature of loss pro- 
cesses. The fall in luminosity is exponential, on a character- 
istic ~ 10* yr timescale. For our fiducial Lmhd ~ 100 kpc 
simulation, I/i.4ghz falls by an order of magnitude in sev- 
eral hundred Myr; the decrease is faster for smaller values of 
Lmuo, which corresponds to strong damping. If only non- 
linear Landau damping operates, the decline in luminosity 
is very slow, and insufficient to turn off radio halos. We also 
show how Li.4GHz evolves if we ignore diffusion (i.e., the no 



damping limit) or adiabatic losses in equation (38 1. 

The streaming speeds relative to the wave frame for 
100 GeV CRps at a radius of 100 kpc are shown in fig- 
ure [3};, for different values of Lmhd and if only no n- linear 
Landau damping dominates. We see that even if the stream- 
ing speeds start out slow, they can quickly become super- 
Alfvenic as the CR density drops. This non-linear behavior, 
which is due to the unusual k oc 1 /V/ scaling of the diffusion 
coefficient for turbulent damping, allows very fast streaming. 
It is not seen if only non-linear Landau damping operates; in 
that case, hd ~ 0{v\) at all times. Note from equation ([7| 
that I'D — t'A ~ (A/3I/z)c. Thus, as hd — >■ c, A — >■ Lz, and our 
equations break down, as the CRs can no longer be described 
by a distribution function. Instead, a fully kinetic approach 
is needed. This limitation is relatively unimportant since by 
this stage CRs are no longer self-confined but stream freely 
along field lines; thus, turn-off is extremely rapid. 

The individual contributions to /p at 100 GeV are 
shown in Figure [3|i in the Lmhd ~ 100 kpc case. The val- 
ues are taken at a fixed radius of 100 kpc, and displayed as 
a function of time. Interestingly, no one process dominates 
(and we have verified that adiabatic and diffusive losses act- 
ing in tandem are much more effective than either process 
alone). Initially, adiabatic losses dominate, although they 
decrease continuously with time. This is to be expected, 
since the adiabatic loss term is proportional to /p, and de- 
creases as /p falls. On the other hand, the diffusion loss 
term from turbulent damping is independent of /p (equa- 
tion (28l), and thus independent of time as long as there is 
a spatial gradient. At f ~ 190 Myr, the profile at r ~ 100 kpc 
flattens (see Fig. [3^), and all terms plummet, although the 
adiabatic loss term falls most drastically. In i|5] we explore 
the nature of this change when the profile flattens: inside the 
flat core, /p changes since it is determined solely by the flux 
at the outer boundary of the core. From this plot, we can 
see why adiabatic and diffusive losses in tandem are much 
more efficient that either alone: adiabatic losses are much 
more effective in the early stages when the profile is cen- 
trally peaked, while diffusive losses are more effective once 
the profile fiattens. We also see that inward turbulent advec- 
tion is non-negligible but subdominant. It also changes once 
the region becomes incorporated inside the flat core (since 
once again only the flux through the core boundary matters 
at that point). 

Fig[3]3 shows 47rp^fp (i.e., CR density) for 100 GeV CRs. 
As previously discussed, the CR density profile develops a 
flat core, which expands in size at roughly the streaming 



^ This is mostly due to projection effects; note that the CR radial 
profile does evolve significantly; see Fig. |3^. 



speed. Meanwhile, the normalization of the CR profile falls 
continuously, even for the flat portion. The end result is 
a proflle in which the CR proflle has completely flattened 
and fallen by several orders of magnitude by the end of the 
simulation. 

Finally, the expected 7-ray flux as calculated from equa- 
tion ( 65 1 is shown in Fig [sj. Observed upper limits are 
also shown; note that our initial conditions are consistent 
with these upper limits. Due to the flnite momentum grid 
Pp ^ 5000, our calculations are only accurate in the range 
E^i < 200 GeV, although the high energy CRs stream so 
quickly that the CR transport equation quickly breaks down, 
in any case. The gamma-ray fluxes decline extremely rapidly 
with time, with the decline being much sharper at higher en- 
ergies, due to the fact that higher energy CRs stream faster. 
The upshot is that at the ~ 0.3 - 1 TeV [Ecn ~ 3 - 10 
TeV) energies probed by imaging air Cherenkov telescopes 
(MAGIC, HESS, VERITAS), the decline in gamma-ray lu- 
minosity is very rapid. Any detection of gamma-ray emis- 
sion at these energies, where a source is not immediately 
apparent (suggesting that it is long-lived), would strongly 
disfavor the model of CR streaming presented here. How- 
ever, the £;.^ ~ 0.1 - 3 GeV {Eck ~ 1 - 30 GeV) energies 
probed by Fermi correspond to CRs which stream and turn 
off gamma-ray emission more slowly. The latter is thus a 
more robust measure of the cluster's CR injection history. 
Note that since {E-y) ~ 1/8 (-Ecr), gamma-ray emission at 
E^ ~ 10 GeV corresponds to the Eck ~ 100 GeV CRs rel- 
evant for ~GHz radio emission, and declines by a similar 
amount . 

By the same token, the energy dependence of CR 
streaming implies that radio luminosity turns off more 
slowly at lower frequencies. We show this in Fig[7| We can 
also see this in figure |4] which plots the distribution function 
versus momentum. The higher energy CRps drop in den- 
sity much faster than lower energy CRps. The correspond- 
ing high energy synchrotron emission then also drops faster. 
This behavior could explain radio halos such as Abell 521, 
which is detected at 240, 325 and 610 MHz, but not at 1.4 
GHz, implying a cutoff or strong spectral curvature at high 
frequencies ( Brunetti et al.|2008 |. We therefore predict that 
at the low frequencies probed by LOFAR, radio halos should 
be significantly more abundant. 



4.2 Coma Cluster 

We now turn to Coma, a prototypical giant radio halo. We 
focus on the differences with our previous example, Perseus. 
Due to the flat and extended observed surface brightness 
profile, which extends out to 1 Mpc (Fig |5^, observations 
from [Deiss et al.| ( |1997[ )), the inferred CRp distribution 
is much flatter. Indeed, for the B-field we have assumed, 
B oc p"^, Ob ~ 0.3 (which is consistent with rotation mea- 
sure observations (iBonafede et al. 



, the radio proflle 

can be fit by a nearly flat CRp density (Fig. [5^). This large 
flat inferred profile is suggestive that extensive streaming 
has already taken place. Coma thus presents an interesting 
challenge, to see if a significant decline in luminosity is pos- 
sible despite the absence of significant CR gradients except 
a small one at the outer boundary. We emphasize once again 
that for turbulent wave damping, our solutions are indepen- 
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Figure 3. Simulation results for the Perseus cluster, (a) Radio surface brightness of Perseus for Lmhd = 100 kpc. Observations from 
[Pedlar et al.| | [T990[ |. (b) The time evolution of the Perseus cluster's radio luminosity for different levels of damping. The solid lines show 
MHD turbulence damping at various strengths. The dashed line shows non-linear Landau damping, (c) Cosmic ray streaming speeds of 
ICQ GeV CRs at a fixed radius of 100 kpc. (d) Different contributions to /p for the Lmhd = 100 kpc Perseus simulation, (e) Radial 
distribution of 100 GeV protons for Lmhd = 100 kpc. (f) Predicted gamma-ray fluxes. Upper limits from observations are at higher 
energies than those plotted here. 
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Figure 4. CR distribution function versus momentum at a fixed 
radius of 100 kpc. The dropoff time scales with energy, leading to 
the spectral steepening discussed above. 



dent of the magnitude of the CR gradient V/. Our solutions 
depend only on where V/ is non-zero, and its sign. 

We show our results in Fig. [5] All figures are analogs of 
those for Perseus in Fig. [s] (except we adopt r=300 kpc as our 
fiducial radius when displaying time-varying quantities — 
due to the much larger extent of the Coma radio halo, 
this is a more representative radius), and we again adopt 
Lmud = 100 kpc for our fiducial model. In Fig [5^, we see 
that the surface brightness falls in normalization, but does 
not significantly change shape, as for Perseus. The decline 
in I/1.4GHZ is slower than for Perseus, although Z/i.4GHz is 
still down by an order of magnitude after ~ 600 Myr for the 
fiducial model, and decHnes more quickly with more vigorous 
turbulence as expected. Non-linear Landau damping alone 
produces a slow decline. While streaming is super-Alfvenic 
(Figjs];), it does not 'run away' with time as quickly as for 
Perseus. The acceleration of CR streaming is tied to the 
decline of the CR density, which is slower in this case. 

Since the flat region encompasses the entire cluster at 
the outset, there is no transition in energy loss regimes as for 
Perseus (where the proflle gradually flattens). Instead, loss 
rates vary mildly with time (Fig[5|i), with diffusive losses al- 
ways more important than adiabatic losses, which are essen- 
tially negligible. This can be understood from the fact that 
adiabatic flux at the outer boundary scales with /p(-Rmax), 
which is small (see ijsjfor more discussion), whereas the dif- 
fusive flux is independent of /p(i?max). The flat CR density 
profile simply decreases in normalization with time (Fig[5]5). 
Similar to Perseus, Coma's gamma-ray flux declines quickly 
with time, particularly at the high energies associated with 
imaging air Cherenkov telescopes. 

Spectral steepening in Coma's radio emission has been 
seen in multi- frequency observations ( Brunetti et al.|[2012 |, 
a feature which occurs naturally in our models due to the 
energy dependence of CR streaming. We show this in Fig 
[6] spectral steepening very similar to that observed arises. 
Given the flat inferred proflle of Coma, which suggests that 
substantial streaming has already taken place, this raises 



the possibility that a power-law population with a slightly 
higher normalization was transformed by streaming into the 
curved population we see today. 

We have chosen a rather extreme case of a completely 
flat proflle, to illustrate that radio halo turn-off is still possi- 
ble in this case. The observational data also permit an initial 
CR proflle that is less flat. Using the same B-fleld, we can 
still reproduce the observations very well with a proflle that 
has a mild central peak. Since we now have a signiflcant den- 
sity gradient, the radio luminosity can drop off faster in the 
beginning, although the overall evolution is qualitatively the 
same as before. The streaming speeds ramp up faster than 
in the flat profile fit. 

We have assumed a magnetic fleld proflle B oc p"^, 
with qb = 0.3. This choice assumes the B-field is in rough 
equipartition with turbulence, as discussed in §2.4| agrees 
with Faraday rotation measures, and enables us to repro- 
duce the observed surface brightness distribution. However, 
the Faraday rotation measurements are consistent with a 
range of values q:b ~ 0-4 — 0.7 at la ( Bonafede et al.|2010 |. 
A steeper scaling of B with density implies lower B-fields at 
the cluster outskirts; in this case the rate of turn-off and de- 
velopment of spectral steepening will be slower. While low- 
ering B decreases the Alfven speed, the dominant effect is a 
decrease in the CR flux F due to streaming, which scales as 
^3/2y£^i/^^ (equation (28 1). The rate of CR streaming is set 



by the minimum value of this flux, which generally occurs at 
the cluster outskirts. As we shall we shall see in ^ in this 
regime we can approximate /p ~ 3F{Rf ,p) / Rf . Specializ- 



ing to our model for Coma {Anp fp{t — 0) 



IQ- 



for p = 100, nc(0)/nc(l Mpc) 



17, and Bo = 5 ^G): 



to 



370 Myr L 



1/2 



^70 



ciB-O.S 



(77) 



1 Mpc 



Thus, a steeper scaling as = 0.7 would not permit turn 
off on an acceptably short timescale for a very ffat profile (it 
could still be possible for a less flat proflle, as above, but note 
that the observations cannot be flt well by a centrally peaked 
CR profile with this steep B-fleld scaling). Note that radio 
relic measurements are consistent with strong, ~ /xG fields 
at the cluster outskirts ( Feretti et al.|2012 l; in addition, the 
very strong turbulence at the cluster outskirts could be con- 
sistent with lower values of Lmhd than the constant value 
we have assumed. As we reiterate in the Conclusions, com- 
plex issues regarding magnetic fleld topology and strength 
are best further explored with 3D MHD simulations. 



5 ANALYTIC EXPRESSIONS 

In certain limiting cases, the evolution of the CR popula- 
tion can be derived analytically. These solutions serve two 
purposes: they serve as tests of our numerical code, partic- 



ularly the regularization scheme (i3.1|, and they also give 



physical insight into the behaviour of our solutions, and the 
circumstances under which particular processes dominate. 

In the absence of sources Q and ignoring the negligible 
Coulomb and hadronic losses, we can write the CR transport 



equation ( 38 1 as 



-V-F 



Dip _a/p 

Dt ^ dt 

where the total CR flux F = Fadia + Fs 



(78) 

+ Fturb, is made 
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Figure 5. Simulation results for the Coma cluster, (a) Radio surface brightness of Coma for Lmhd = 100 kpc. Observations from |Deiss| 
|et al.| ( [1997^ . (b) The time evolution of the Coma cluster's radio luminosity for different levels of damping. The solid lines show MHD 
turbulence damping at various strengths. The dashed line shows non-linear Landau damping, (c) Cosmic ray streaming speeds of 100 
GeV CRs at a fixed radius of 300 kpc. (d) Different contributions to /p for the Lmhd = 100 kpc simulation, (e) Radial distribution of 
100 GeV protons for Lmhd = 100 kpc. (f) Predicted gamma-ray fluxes. Upper limits are taken from | |Arlen et al.|2012l l with a = 2.5. 
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Figure 6. Luminosity as a function of energy for the Coma sim- 
ulation including observations from | |Brunetti et al.||2012^ . The 
momentum dependence of the streaming speed leads to a spec- 
tral steepening very similar to observation. 
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Figure 7. Luminosity dropoff in Coma for different frequencies. 
High energy CRs stream more quickly, so the higher frequency 
signals drop faster. 



up of the fluxes due to adiabatic losses in the wave frame, 
streaming relative to the wave frame, and turbulent advec- 
tion respectively. We have approximated the Lagrangian 
derivative by the Eulerian derivative D/p/Dt « dfp/dt, 
since va • V/p is initially small and becomes increasingly 
negligible as the profile flattens. 

As we have seen, the CR profile generally develops a 
flat inner core within some radius Rj , outside of which it 
declines. The flat core stems from the fact that while V ■ F 



increases inwarcf"^ an inverted CR proflle cannot develop, 
since CRs cannot stream up a gradient. Thus, a flat core 
develops, while its normalization and radius Rf evolves due 
to the net flux of CRs from its outer boundary. In particular, 
if we set F = Fr and integrate equation ( 78 1 over the volume 
of the flat region, we obtain: 

^7vRj,URf,p) = -A7vR^F{Rf,p) 



fARf:V) = 



SF{Rf,p) 
R 



(79) 



where we have used the fact that fp{r,p,t) is independent 
of r for r < Rf, and the divergence theorem. The evolution 
of the entire profile can then be described by 



fp{r,p,t) = 



3F{Rf,p,t) 
Rf(p,t) ' 

-V-F(r,p,f), 



r < Rf{p,t) 
r > Rf{p,t) 



(80) 



where the "flatness front" Rf{p,t) is determined from 
fp{0,p,t) = fp{Rf,p,t), or: 

f /n m r 3F{R({p,t'),p,t') 

- / \/-F{Rt{p,t'),p,t')dt' (81) 
Jo 

As we have seen, Fstroam and Fadia are the most impor- 
tant fluxes, while Ftmb is subdominant. Let us now consider 
the limiting cases when only one is at play. 

Cosmic-Ray streaming only. We have: 



Fgtri 



(82) 



in^p^mfloVA 

Since -Fstream is independent of fp and depends only on 
plasma parameters (specifically, the B-field, turbulence and 
density profiles), in our model where the gas properties are 
time-steady (and thus in hydrostatic and thermal equilib- 
rium), Fstraa.in{r,p) is independent of time. Thus, A{r,p) = 
V ■ Fstrcam('', p) IS also tlme-iudependeut, and we have for 
r > Rf{p,ty. 



f{r,P,t) = ,f(r,p,0) - A{r,p)t- r > Rf{p,t) 



(83) 



i.e., the distribution function outside the fiatness front falls 
linearly with time. More generally, we can solve for the fiat- 
ness front Rf and the overall solution both inside and outside 



R{ via equations (801 and (811 



To compare this analytic solution with our simulation 
we ran a simulation for Perseus and for Coma where only 



the diffusion term was used in (381, and all other terms ig- 



nored. The resulting CR densities for Perseus can be seen 
in figure [5] In this plot we show the CR density versus time 
at 100 GeV at a few select radii. The solid curves are the 
simulation, and the dotted lines are the analytic solution for 
a non-fiat profile. The match is essentially perfect - the den- 



sities decrease at a constant rate (equation (83l) until the 



fiatness front catches up to each radius. After this point the 
densities follow the same single curve corresponding to the 



This condition holds as long as F increases more slowly than r. 



In our case, the dominant fluxes _Fatr 



^3/2 



oc p' 



,3(13/2 



(equa- 



tion ( |82[ l) clearly increases inward, and -Fadia fv\ is at most 
fiat or increases inward. Thus, |V ■ F| clearly increases inward. 
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from CR streaming Fstream (equation | |82| ) is important. The 
dotted lines (bold curve) show the analytic solution for outside 
(inside) the flat front respectively; the solution initially follows the 
dotted curves until it intersects the green curve, when it follows 
the flat front solution. The solid lines show the simulation results, 
which match the analytic solution almost perfectly. 
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Figure 9. The same as for Fig[8]but for Coma. Since the profile 
is already almost flat, this is a test of the /p = —SFcriRi) / Ri = 
constant regime. 



evolution of the flat region. The same results for Coma show 
the agreement in the regime when the profile is already flat. 
In figure [9] the CR profile in Coma is nearly flat to begin 
with. Before very long the profile is flat across the entire 
simulated space and the fp = —3Fci{Ri)/Ri regime kicks 
in. Again, the agreement between simulation and analytic 
solution is very good. This implies that our regularization 
of the CR streaming term (which is needed to prevent un- 
physical oscillations with such flat profiles) is not so strong 
that it artificially changes the rate of diffusion. 
Adiabatic expansion only. We have: 



1 df 
3 op 



(84) 
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Figure 10. CR densities for Perseus in the absence of any diffu- 
sion, i.e. only adiabatic losses are used. The dotted line represents 
the analytic solution i85\ and should be compared to the blue line. 



tribution function fp and hence is time-dependent. We can 
readily solve this in the approximation that va oc p°B-0 5 ^ 
const (since it varies very weakly with radius), and fp oc 
p^" , approximately independent of radius. Then, for r ^ _Rf, 
we have f = V ■ F^dia ~ ~2avAfp/3r, or: 



fp{r,p,t) ^ /p(r,p,0)exp 



2avAt\ 



(85) 



Thus, outside the fiatness front, the distribution function 
falls exponentially with time, with e-folding time of order 
the Alfven crossing time (which becomes long at large radii). 
To solve for the evolution of the fiatness front and the entire 



profile, we insert equation ( 85 I into equation ( 84 1 and hence 



equation (801 and (811. Note that we are only required to 
evaluate the fiux F^dia for r ^ Rf, where equation ( |85[ ) is 
valid. We compare this analytic expression with a Perseus 
simulation that has no diffusion in Fig |10[ The dashed line 
depicts (|85| for r = 100 kpc and p = 100. Ahhough the fit 
isn't perfect, the simulated values do fall exponentially with 
time until the flatness front catches up, with e-folding time 
comparable to that determined from l\85\ . This is perhaps to 



Thus, unlike the preceding case, the flux depends on the dis- 



be expected, since ( 85 I assumes that the quantity afp does 
not vary significantly with radius, which is not typically the 
case. 

These solutions allow us to understand the nature of the 
numerical solutions we previously obtained. Coma, where 
the initial profile is almost completely flat, is obviously in 
the fp = -3F{R^!,^)/Rf regime; moreover, Fstream(-R/) > 
Fadia(^?f ), since the latter scales with the (small) value of the 
distribution function at the outer boundary. The evolution 
of the flatness front in Perseus is more interesting. Initially, 
even though vo — va ^ 0{va), adiabatic losses dominate, 
since the distribution function falls exponentially with time 
(rather than linearly with time, for streaming losses). How- 
ever, Fadia oc / also falls exponentially with time, while 
Fstream is Independent of time. Thus, streaming losses will 
always dominate at late times. Equivalently, the velocity as- 
sociated with adiabatic losses, va, is constant with time. 
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while the streaming velocity vd oc l//p increases with time: 
as the number density of cosmic rays fall, the confining wave 
amplitude SB/B falls, and cosmic rays can stream progres- 
sively faster. 



6 CONCLUSIONS 

Shocks generated during hierarchical structure formation 
are expected to accelerate cosmic rays via diffusive shock 
acceleration. These cosmic rays in turn interact hadroni- 
cally with thermal nucleons to produce pions, which de- 
cay to produce relativistic electrons. Tracking these well- 
understood processes, and assuming magnetic fields given 
by Faraday rotation measurements, leads to predictions for 
radio halo emission consistent with those observed ( |Pfrom-| 
mer||2008 1. However, this model predicts that every cluster 
hosts a bright radio halo. This is at odds with the observed 
bimodality of cluster radio emission: the majority of clus- 
ters are radio-quiet, and an order of magnitude fainter than 



the radio-loud population (Feretti et al. 2012 Brown et al 



|2011[ ). Radio loudness is strongly associated with merger ac- 
tivity. For this reason, the turbulent re-acceleration model 



(Brunetti et al. 2001 Petrosian 20011, where this assocl 



ation occurs naturally, is often favored. However, this still 
begs the question as to why hadronically induced radio emis- 
sion is not omnipresent. All of the associated physics is well 
understood, and at face value the observations then require 
that CRp acceleration efficiencies be reduced by an order of 
magnitude below canonical valuej^ 

(20111 took an important step forward 



Enfilin et al. 



when they suggested that CRp's could potentially stream 
super-Alfvenically, turning off radio halos. However, they 
assumed streaming speeds of order the sound speed vo ~ Cs 
instead of calculating il p^ and posited steady-state CR pro- 
files that represent equilibria between outward streaming 
and inward turbulent advection, despite the long timescales 
for equilibration. In this paper, we attempt to place CR 
streaming in clusters on a more rigorous footing, by calcu- 
lating the microphysical streaming speed as a function of 



plasma parameters in the self-confinement picture (Lerche 
19671 iKulsrud & Pearce'1969; Wcntzel 1969, Skilhng 1971| 



We then solve the time-dependent CR transport equation 
(albeit in ID) to see how the radio luminosity evolves with 
time. Our conclusions are as follows: 

• CR streaming speeds depend on the source of wave 



damping. Non-linear Landau damping (e.g., (Felice & Kul 
|srud''2001 



as assumed in 'Enfilin et al. (2011 



l) is too weak 

to sufficiently inhibit wave growth, and vo ~ wa- However, if 
waves are instead damped by turbulent shear ( Yan & Lazar- 
ian|2002' Farmer fc Goldreich|2004 1 , they can be sufficiently 



suppressed that super-Alfvenic streaming «d S> va is possi- 
ble. Moreover, vy) — va (x 7/ncR(> 7), (where 7 is the CR 
Lorentz factor) so that: i) higher energy cosmic rays stream 



more rapidly; ii) CR streaming speeds continually increase 
as ncR declines due to streaming. 

• Streaming relative to the Alfven wave frame can be in- 
corporated into the CR transport equation via a diffusion 
term. For turbulent wave damping, the diffusion coefficient 
K oc 1/V/p (where /p is the distribution function), so that 
remarkably V- («;V/p) is independent of V/p. Thus, CRs can 
continue to stream unabated in giant radio halos (such as 
Coma) despite their fairly fiat inferred CR profiles. Stream- 
ing is still sensitive to the sign of V/p (since CRs can only 
stream down a gradient), and for fiat profiles we must imple- 



ment numerical regularization ( Sharma et al. 2009 1 to ensure 
stable solutions. We test our solver for the CR distribution 
function against a code where CR mediated AGN heating 
is solved in the fluid approximation ( Guo fc Oh||2008" l. The 
solutions are identical. Note that CR heating is unaffected 
by super-Alfvenic streaming, since it scales as va ■ VPc and 
Pc is dominated by ~GeV CRs, where streaming is Alfvenic. 

• CR transport is thus clearly modifled by ICM turbu- 
lence. Besides its effects on wave damping, turbulence can 
also advect CRs so that they roughly trace the gas den- 
sity profile, creating a centrally peaked CR distribution. 
For the mildly subsonic ~ ua ~ 100 km s"'^ turbulence 
we assume, outward streaming dominates inward advection. 
Moreover, this trend increases with the amplitude of turbu- 
lence. It is therefore consistent with the fiat inferred CR pro- 
files in non cool-core clusters, which have generally stronger 
turbulent motions. Such a trend is hard to understand in 



scenarios where turbulence only draws CRs inward (Enfilin 
et al.|2011||Zandanel, Pfrommer fc Prada|20l"2l ) 



We then perform numerical time-dependent calcula- 
tions of CR streaming, assuming an initial profile consistent 
with radio observations at 1.4 GHz. We find that the radio 
luminosity falls by an order of magnitude in several hundred 
Myr, both in a prototypical radio mini-halo (Perseus) and 
giant radio halo (Coma). The latter effect is particularly in- 
teresting in light of the fiat inferred CR profile, and arises 
only for turbulent damping of MHD waves; if only non-linear 
Landau damping is at play, the turn-off is slow. Indeed, the 
inferred fiatness of the CR profile suggests that streaming 
has already been at play in these systems. We also build an 
analytic model which aids in physical understanding. Adia- 
batic losses dominate until the profile fiattens, when diffusive 
losses dominate. The turn-off timescale in the later stage is 
set by the lowest value of the CR fiux F oc B^^^ / Lmhd, 
generally at the cluster outskirts. The energy-dependence of 
CR streaming means that spectral curvature develops, and 
radio halos turn off more slowly at low frequencies, both 



consistent with observations (Brunetti et al. 2008 20121 



Streaming also rapidly diminishes the 7-ray luminosities 
at the E-y ~ 0.3 — 1 TeV energies probed by imaging air 
Cerenkov telescopes (MAGIC, HESS, VERITAS), but not 
for the lower energies ~ 0.1 — 3 GeV probed by Fermi. 
The latter is therefore a more robust probe of the CR injec- 
tion history. 



A bimodality in cluster B-fields, with larger values during 
the turbulent, radio-loud state, appears inconsistent with clus- 
ter rotation- measure observations ( [Bonafede et alT| | |2011| , and 
references therein). 

^■^ In fact, given their assumptions, we find that cosmic rays 
should only stream Alfvcnically. 



The primary contribution of this paper is a physical 
proof of principle for turning off hadronically induced emis- 
sion. Our ID streaming calculations by nature omit impor- 
tant details best clarified by 3D MHD simulations. Chief 
amongst these are the effects of magnetic topology. We have 
effectively assumed radial magnetic fields in our ID calcu- 
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lations. Of course, magnetic topology greatly influences the 
true value of macroscopic transport coefficients. There is 



some evidence both from observations (Pfrommer & Dursi 



20101 and cosmological MHD simulations (Ruszkowski et al 
2011 1 that outside the core, magnetic fields are largely ra- 
dial, driven either by cosmological infall, or the magneto- 
thermal instabil ity (MTl; Balbusl ( |2000| ; |Parrish, Stone~fcl 



Lemaster (20081). Alternatively, turbulence could fully tan- 



gle magnetic fields 


Ruszkowski & Oh||2010 2011 Parrish, 


Quataert & Sharma 


2009 1. CRs have to follow the same field 



lines that thermal particles do, albeit with a larger gyro 
radiuj^ As long as cross-field diffusivity remains small, 
transport coefficients should scale similarly; in the limit of 
a fully tangled field with a coherence length significantly 
larger than the gyro radius, a random walk in 3D rather 
than ID will reduce the diffusion coefficient — >■ Kp/3, just 
as it reduces the Spitzer-Braginskii value for thermal con- 
ductivity by a factor of 3. This will effectively increase all 
quoted timescales in this paper for pure streaming by a fac- 
tor of ~ 3. In the future, it would be interesting to conduct 
fully self-consistent 3D MHD simulations which include CR 
streaming, motivated and guided by the estimates here. We 
have also incorporated the advective effects of gas motions 
only in the diffusive approximation. Coherent bulk motions 
due to mergers or sloshing could potentially have stronger 
effects. More light on the nature of ICM motions in radio- 
bright halos from Astro-H (e.g., [Zhuravleva et aT] ( |2012[ ); 
Shang & Oh ( 2012b|a l) will surely help. We are also agnos- 
tic as to the cause of radio halo turn-on, which is clearly 
related to gas motions stimulated by mergers, and could be 



due to turbulent reacceleration of seed CRe ( Brunetti et al 



2001 Petrosian 20011, inward advection of CRs from the 



cluster outskirts ( [Enfilin et al.|["2011[ ), or perhaps have sep- 



arate mechanisms for different classes of radio halos ( Zan- 



danel, Pfrommer fc Prada,,2012j ). Such issues await clarifi- 
cation from low frequency radio observations by LOFAR. 
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